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PREFACE 


11^ 

y.  / 

The  initial  co-sponsored  Air  Force  Systems  Command/Naval 
Material  Command  Science  and  Engineering  Symposium  was  held  at  the 
Naval  Amphibious  Base,  Coronado  on  16  -  19  October  1978.  The  theme 
of  the  1978  Symposium  was  "Advanced  Technologies  -  Key  to  Capabilities 
at  Affordable  Cost." 

The  objectives  of  this  first  joint  Navy/Air  Force  Science  and 
Engineering  Symposium  were  to: 

.  Provide  a  forum  for  military  and  civilian  laboratory 
scientific  and  technical  researchers  to  demonstrate 
the  spectrum  and  nature  of  1978  achievements  by  their 
services  in  the  areas  of 

.  Armament  .  Human  Resources 

.  Avionics  .  Materials 

.  Basic  Research  .  Propulsion 

..  Flight  Dynamics 

.  Recognize  outstanding  technical  achievement  in  each 
of  these  areas  and  select  the  outstanding  technical 
paper  within  the  Navy  and  the  Air  Force  for  1978 

.  Assist  in  placing  the  future  Air  Research  and 

Development  of  both  services  in  correct  perspective 
and  to  promote  the  exchange  of  ideas  between  the  Navy 
and  Air  Force  Laboratories 

.  Stress  the  need  for  imagination,  vision  and  overall 
excellence  within  the  technology  community,  assuring 
that  the  air  systems  of  the  future  will  not  only  be 
effective  but  affordable. 

Based  upon  the  success  of  the  initial  joint  symposium  (which 
was  heretofore  an  Air  Force  event),  future  symposia  are  planned  with 
joint  Navy/Air  Force  participation. 
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ABSTRACT 


The  Airborne  Electronic  Terrain  Map  System  (AETMS)  Is  a  compact, 
low  cost  system  capable  of  generating  real-time  displays.  It  couples 
directly  to  the  aircraft *s  navigational  computer  and,  In  effect,  auto¬ 
mates  the  map  reading  process.  It  uses  state-of-the-art  components  and 
Defense  Mapping  Agency  data  bases.  The  system  architecture  Is  quite 
general  and  can  be  adapted  for  different  applications;  Its  core  elements 
can  be  easily  shared  In  different  user  configurations. 

The  United  States  Air  Force  Avionics  Laboratory  (AFAL)  has  a  pro¬ 
gram  to  develop  a  brassboard  prototype  of  this  system.  The  feasibility 
of  the  AETMS  hks  already  been  partially  established  by  simulation  of  the 
architecture  on  a  general  purpose  computer  at  AFAL.  This  simulation  has 
been  extended  to  display  terrain  Information  In  near-real- time  on  a  CRT 
linked  to  the  computer.  The  displays  are  generated  from  stored  poly¬ 
nomial  terrain  models  developed  at  the  United  States  Army  Engineer  Topo¬ 
graphic  Laboratory. 

The  flexibility  of  electronic  terrain  maps,  compared  to  paper  and 
film  maps,  lends  Itself  to  many  Air  Force  applications.  Low-level  ter¬ 
rain  displays  are  the  most  obvious,  ^lle  others  Include  prediction  of 
terrain  problems  beyond  the  Immediate  horizon,  new  methods  of  obtaining 
navigational  fixes  at  low  altitudes,  and  terrain  navigational  systems 
not  constrained  to  fixed  flight  profiles. 
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Introduction 


The  USAF  Avionics  Laboratory  (AFAL)  is  currently  developing  an  Air¬ 
borne  Electronic  Terrain  Map  System  (AETMS) .  This  effort  was  originally 
motivated  by  the  desire  to  reduce  map  reading  chores  during  low  level 
missions  in  single  seat  aircraft  (Figure  1) .  In  this  highly  d3mamlc  en¬ 
vironment,  the  task  of  correlating  navigational  data  with  paper  maps  in 
order  to  visualize  terrain  is  more  difficult  than  usual.  This  paper 
briefly  describes  the  approaches,  results,  and  status  of  our  efforts 
toward  automating  the  map  reading  process  for  pilots.  (Reference  1 
gives  a  more  comprehensive  discussion  of  AETMS.) 

Our  primary  goal  is  to  develop  a  black  box  which  will  contain  a 
large  area  map  along  with  Specialized  hardware  and  software  for  acces¬ 
sing  the  terrain  data  in  real-time.  Ground-based  computers  can  fulfill 
this  function,  but  no  such  system  exists  with  the  size,  speed,  and 
flexibility  required  for  economical  employment  in  the  airborne  environ¬ 
ment  . 


Our  second  goal  is  to  develop  an  airborne  display  system  which  will 
generate  dynamic  terrain  displays  using  output  data  from  the  electronic 
map.  Even  though  the  electronic  map  can  be  used  to  generate  a  great 
variety  of  displays,  the  airborne  display  generator  must,  in  order  to 
keep  complexity  at  a  minimum,  use  simple  algorithms  to  address  the  prob¬ 
lem  of  real-time  conversion  of  streams  of  elevation  data  into  display 
signals.  The  results  should  not  be  expected  to  achieve  state  of  the  art 
realism  in  computer  generated  Imagery;  however,  computer  simulation 
tests  have  produced  some  remarkably  good  terrain  displays  with  the  fast, 
simple  algorithms  suitable  for  airborne  systems. 

Terrain  Modeling 

The  standard  level  1  digitized  terrain  data  base  developed  by  the 
Defense  Mapping  Agency  (DMA)  Is  available  and  adequate  for  an  airborne 
map  system.  The  standard  format  is  600  elevation  values,  each  com¬ 
prised  of  16  bits,  to  represent  approximately  one  square  mile  (9600  bits 
per  square  mile) .  It  is  readily  apparent  that  large  areas  require  huge 
amounts  of  data.  Attempts  to  reduce  this  amount  by  thinning  out  the  ma¬ 
trix  of  elevation  values  result  In  Increased  errors. 

As  an  alternate  approach,  elevation  values  are  not  stored  directly; 
Instead,  locally  valid  analytic  surface  models,  polynomials,  are  used. 
The  polynomial  modeling  techniques  developed  at  the  Army  Engineer  Topo¬ 
graphic  Laboratories  (ETL)  use  smoothly  joining  surface  functions  which 
reduce  the  DMA  data  without  reducing  data  accuracy  (Reference  2) .  In 
this  approach,  polynomial  co-efficients  are  stored  in  a  digital  memory 
and  these  can  then  be  used  to  compute  the  elevation  of  any  point. 
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Imagine  a  two-dimensional  coordinate  grid  covering  a  map  surface 
and  sectioning  the  surface  Into  patches  of  approximately  0.2  by  0.2 
miles.  In  other  words,  each  square  mile  Is  covered  by  a  set  of  25 
patches.  The  corners  of  these  patches  (or  the  Intersections  of  grid 
lines)  are  called  nodes.  In  the  polynomial  models,  a  set  of  co-effl- 
clents  Is  associated  with  each  node  and  characterizes  a  bivariate  poly¬ 
nomial  centered  at  the  node.  This  bivariate  polynomial  has  a  domain  of 
validity  over  the  four  adjacent  patches  having  the  node  as  a  common  cor¬ 
ner. 

Given  a  desired  sample  position  within  a  patch,  the  data  at  Its 
four  nodes  must  be  weighted.  These  weights  ate  functions  of  the  rela¬ 
tive  position  (patch  coordinates)  of  the  sample  with  respect  to  each  of 
the  four  nodes.  The  advantage  In  data  reduction  of  using  polynomials 
Is  paid  for  by  an  Increase  In  computational  complexity.  For  discrete 
values,  the  weighted  average  Is  equivalent  to  a  bilinear  Interpolation. 
With  pol3momlals,  each  of  the  four  nodal  polynomials  must  be  evaluated 
at  the  sample  point  and  then  averaged  to  yield  a  valid  elevation  value. 

Current  working  estimates  for  the  airborne  data  base  are: 

25  patches  per  square  mile. 

34  bits  per  set  of  four  polynomial  to-6fflclents 
(a,  b,  c,  d),  where  the  bivariate  polynomials  have 
the  form  a  +  bx  +  cy  +  dxy. 

6  bits  per  patch  to  encode  cultural  features. 

This  Is  equivalent  to  1,000  bits  per  square  mile.  Including  cultural 
data.  This  Is  large,  but  significantly  smaller  than  the  standard  den¬ 
sity  of  9,600  bits  per  square  mile  for  encoding  elevation  data  only. 
On-going  studies  at  ETL  suggest  that  even  further  reductions  may  be  pos¬ 
sible.  DMA.,  ETL,  and  AFAL  are  cooperating  In  further  improvements  to 
the  airborne  data  base. 

Electronic  Map  Implementation 

Memory  technology  provides  a  choice  of  high  density  memory  packages 
In  which  to  store  the  data  base.  Storage  of  a  very  large;  l.e.,  global, 
terrain  map  may  even  require  a  combination  of  memories,  such  as  a  mag¬ 
netic  bubble  memory  backed  up  with  archival  data  on  magnetic  tape.  Bub¬ 
ble  memory  is  particularly  attractive  because  It  is  relatively  inexpen¬ 
sive,  rugged,  and  non-volatile.  Current  16  megabit  prototype  packages, 
developed  for  AFAL  by  Texas  Instruments,  have  storage  capacity  for  a 
regional  map  of  16,384  square  miles.  Magnetic  tape  could  be  used  to 
periodically  update  such  a  bubble  memory.  Programs  are  under  way  to 
significantly  Increase  the  packing  density  of  the  bubble  memory  in  the 
early  1980’ s. 
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The  airborne  system  is  designed  to  achieve  a  sampling  rate  of  one 
million  elevation  samples  per  second.  This  rate  should  be  adequate  for 
most  real-time  applications.  Achievement  of  this  rate  in  a  compact, 
low  cost,  airborne  system  operating  in  a  dynamic  environment,  however, 
presents  several  major  difficulties. 

To  sample  every  microsecond,  four  polynomials  need  to  be  computed 
and  averaged  every  microsecond.  This  is  achievable  only  with  a  cus¬ 
tomized  parallel  processor  known  as  a  polynomial  generator.  The  major 
problem  here  is  the  reconciliation  of  the  slow  access  times  inherent  in 
mass  memories  with  very  fast  parallel  data  transfer  and  processing 
(Figure  2) .  Block  access  time  into  bubble  memories  is  one  or  two  mil¬ 
liseconds,  which  is  reasonably  good  time  for  mass  memories,  but  too 
slow  to  drive  the  polynomial  generator. 

We  can  speedily  access  data,  however,  if  it  is  temporarily  stored 
in  a  dual  port  fast  memory  buffer  with  a  parallel  output  data  channel. 
This  fast  memory  buffer  acts  as  a  dynamic  window  for  the  global  model 
in  mass  memory  and  contains  all  relevant  patch  data  in  the  immediate 
vicinity  of  the  moving  aircraft. 

The  design  of  this  dynamic  buffer  Involves  the  reconciliation  of 
opposing  goals  and  requirements;  namely,  (1)  providing  sufficient  agi¬ 
lity  to  enable  the  buffer  to  adapt  and  respond  to  signals  describing 
curvilinear  trajectories,  (2)  minimizing  I/O  interference  between  the 
updating  and  accessing  processes,  (3)  keeping  the  buffer  capacity 
within  the  minimum  necessary  to  insure  the  containment  of  a  preselected 
usable  window  area.  Knowing  the  characteristics  of  trajectory  dynamics, 
one  can  always  select  a  buffer  memory  of  sufficient  size  to  Insure  ful¬ 
fillment  of  (1)  and  (2) . 

To  avoid  significant  and  undesirable  pauses  (or  down  times)  in  data 
accessing,  some  excess  buffer  capacity  is  always  required  in  a  dynamic 
buffer.  As  an  extreme  measure,  the  buffer  size  can  approach  the  size  of 
the  mass  memory.  The  fast  buffer  memory  system  designed  at  AFAL  repre¬ 
sents  a  significant  breakthrough  for  the  airborne  map  and  display  system 
in  Figure  3.  The  buffer’s  capacity  can  be  minimized  with  respect  to  any 
given  window  size.  In  addition,  it  can  be  economically  constructed  with 
off-the-shelf  components  and  with  ordinary  fabrication  and  packaging 
techniques. 

The  minicomputer  executive  depicted  in  Figure  3  is  generally  res¬ 
ponsible  for  low  frequency  chores,  such  as  external  input  to  the  mass 
memory,  interaction  with  the  navigational  computer,  tracking  the  dynamic 
window,  and  Initializing  the  update  and  scan  controllers.  Once  initia¬ 
lized,  the  update  controller  is  responsible  for  refreshing  the  buffer 
memory  due  to  shifts  in  the  window  position. 
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Each,  scan  pattern  generated  by  the  electronic  map  consists  of  sam¬ 
ples  taken  along  a  set  of  line  segments  which  are  specified  by  the  mini¬ 
computer.  The  scan  controller  determines  map  coordinates  for  each  sam¬ 
ple,  the  buffer  address,  local  patch  coordinates,  and  weights  for  each 
node.  For  each  sample  point.  It  Initiates  transfer  of  buffered  data 
over  a  parallel  data  channel  directly  to  the  parallel  processdr  and  also 
sends  the  patch  coordinates  and  weights  over  separate  data  channels. 

Data  at  four  neighboring  nodes  are  sequentially  accessed  and  processed 
In  order  to  reduce  the  number  of  processing  elements  and  size  of  data 
channels.  (Presently,  additional  circuitry  to  generate  samples  along 
non-llnear  segments  Is  not  justified  by  any  known  requirements.) 

Both  the  update  controller  and  scan  controller  have  special  elec¬ 
tronic  circuits  to  perform  their  respective  functions  automatically, 
which  greatly  eliminates  data  management  problems.  Together,  they  han¬ 
dle  the  high  frequency  calculations  for  the  system.  The  memory  control 
system  has  been  carefully  Integrated  so  that  both  controllers  share  con¬ 
tinual  access  to  the  fast  buf;f er  memory  and  the  scan  controller  can  ac¬ 
cess  all  data  at  a  grid  node  with  a  single  buffer  address. 

The  polynomial  generator  architecture,  also  designed  at  AFAL,  can 
accept  all  data  Into  Its  first  level  of  processing  elements  and  compute 
weighted  averages  In  pipeline  fashion.  A  stream  of  elevation  values  can 
be  directed  to  various  application  modules,  such  as  the  display  genera¬ 
tor  shown  In  Figure  3. 

To  achieve  a  one  megahertz  throughput  rate,  memory  cycle  and  para¬ 
llel  multiply  times  of  250  nanoseconds  are  required.  With  faster  com¬ 
ponents,  throughput  rates  Increase  proportionally. 

Displays 

Figure  4  illustrates  two  typical  sampling  patterns  for  display  pur¬ 
poses.  Sample  spacing,  Ap,  need  not  be  uniform  as  depicted  In  the  dia¬ 
grams.  Position,  orientation  (a  ),  or  type  of  pattern  makes  little 
difference  to  the  data  management^techniques  employed  in  the  memory  con¬ 
trol  system.  In  other  words,  the  time  necessary  for  completing  a  scan¬ 
ning  pattern  Is  proportional  to  the  number  of  samples  In  the  pattern. 

Either  the  radial  or  the  parallel  scan  format  can  be  used  to 
generate  a  moving  contour  display  In  which  brightness  or  color  varies 
with  surface  elevation  (Figure  5) .  This  contour  display  can  be  compared 
with  sensor  data  and  used  to  fill  In  the  radar  shadows.  It  enables  the 
pilot  to  look  ahead,  "see"  behind  hills,  and  plan  his  maneuvers.  Natu¬ 
rally,  one  can  combine  other  Information  with  this  display,  such  as 
navigation  and  target  cueing  Information. 
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The  radial  scan  pattern  has  been  used  to  generate  perspective  dis¬ 
plays,  which  are  more  graphic  than  the  contour  display  and  also  more 
difficult  to  generate  in  real-time.  A  simple  masking  algorithm  can  de¬ 
termine,  by  comparing  successive  samples  along  a  radial,  whether  or  not 
a  sampled  position  is  masked  by  intervening  terrain.  It  can  also  flag 
points  just  on  the  verge  of  being  masked.  These  points  are  used  to  de¬ 
pict  the  ridge  lines  in  the  upper  diagram  in  Figure  6.  This  austere 
format  is  suitable  for  a  heads-up  display.  Ridge  lines  above  the  air¬ 
plane  symbol  are  above  the  aircraft’s  altitude.  In  the  bottom  format, 
sky  shading  and  visible  range  llneis  (contouring)  provide  additional  in¬ 
formation  about  the  surface. 

The  dynamic  contour  and  perspective  formats  are  typical  of  the  dis¬ 
plays  which  have  been  explored  at  AFAL  with  computer  simulations.  Cur¬ 
rently,  software  is  successfully  simulating  some  of  the  buffer  manage¬ 
ment  techniques  intended  for  the  airborne  system  on  a  general  purpose 
PDF  11/45.  This  software  enables  us  to  continually  access  a  DMA  data 
base  stored  on  a  magnetic  disk  and  generate  near  real-time  dynamic  dis¬ 
plays  whose  altitude  and  heading  are  controlled  interactively  with  a  joy 
stick.  With  the  current  system  configuration,  the  PDP  11/45  is  com¬ 
puting  display  symbology  at  a  faster  rate  than  can  be  transmitted  over 
the  data  channel  to  a  RAMTEC  9300  display  refresh  memory.  Plans  are 
under  way  to  Improve  our  throughput  processing  rate  with  a  specialized 
refresh  memory  system  directly  coupled  to  a  peripheral  MAP  200  Array 
Processor. 

The  airborne  display  generator  will  Implement  in  hardware  the  real¬ 
time  algorithms  being  developed  at  AFAL.  Its  design  will  be  finalized 
after  interactive  (hands-on)  simulation  experiments,  human  factors  analy¬ 
sis,  and  flight  testing.  The  exploration  of  optimal  display  formats  is 
being  pursued  in  cooperation  with  ETL  and  the  Air  Force  Aerospace  Medical 
Research  Laboratory. 

DMA  is  providing  a  digital  model  of  the  moving  terrain  belt  system 
being  added  to  the  avionics  simulation  complex  at  the  Air  Force  Avionics 
Laboratory’s  System  Avionics  Division.  This  will  provide  a  unlq%e  faci¬ 
lity  for  testing  pilot  reactions  to  various  display  formats  in  a  realis¬ 
tic  cockpit  simulator. 

A  major  program  is  under  way  to  design  and  build  a  brassboard  pro¬ 
totype  of  the  airborne  system  within  the  next  two  or  three  years.  This 
will  be  evaluated  at  AFAL  and  then  undergo  limited  flight  testing. 

Conclusion 

We  have  presented  a  brief  survey  of  a  versatile  and  economical  Air¬ 
borne  Electronic  Terrain  Map  System  whose  high  speed  omnidirectional 
sampling  capability  can  support  many  different  applications.  Some 
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applications  depend  solely  on  navigational  computer  outputs.  Indepen¬ 
dent  of  externally  radiating  sensors  aboard  the  aircraft,  while  others 
can  combine  or  correlate  Information  from  onboard  sensors  such  as  radar 
with  the  terrain  map  Information. 

The  use  of  the  electronic  map  In  two  types  of  terrain  avoidance 
displays,  contour  and  true  perspective,  has  been  described.  The  true 
perspective  display  will  give  a  pilot  an  Immediate  visual  picture  of 
the  ground.  True  perspective  plus  additional  flight  symbology  on  a 
heads-up  display  would  be  a  valuable  aid  In  conditions  of  restricted 
visibility,  especially  In  the  absence  of  radar.  The  system  can  allow 
the  pilot  to  disengage  the  navigational  computer  and  Input  his  choice 
of  state  vectors  Interactively;  l.e.,  fly  the  display.  The  ability  to 
preview  a  dynamic  picture  of  the  terrain,  along  with  target  cueing  sym¬ 
bology,  would  be  useful  In  planning  loif-level  strikes  called  while  In 
the  air. 

Improved  synergistic  systems  combining  the  unique  features  of  the 
electronic  terrain  map  with  various  avionics  sensors  are  foreseeable. 
This  Is  particularly  true  In  low-level  environments  where  much  of  the 
terrain  may  be  above  the  flight  plane.  An  obvious,  but  useful,  example 
Is  the  combination  of  terrain  avoidance  radar  with  the  contour  display 
which  could  provide  the  missing  Information  In  areas  of  radar  shadow. 
Terrain  map  data  plus  radar  altimeter  could  provide  Improved  algorithms 
for  terrain  following  and  contour  navigation  systems  which  are  not  con¬ 
strained  to  pre-established  ground  tracks.  Techniques  for  generating 
true  perspective  could  be  used  to  predict  approximate  radar  Images. 

These  predicted  Images  plus  actual  radar  Images  plus  pattern  recogni¬ 
tion  techniques  may  enable  radar  navigational  fixes  to  be  automated. 

What  Is  applicable  to  piloted  aircraft  also  applies  to  ground-based 
systems.  RPV's  can  be  adversely  affected  by >poor  visibility  or  loss  of 
their  video  communication  channels.  The  electronic  terrain  map  can  be 
used  In  the  ways  described  above  to  enhance  RFV  control  and  guidance 
while  providing  a  ready  back-up  for  the  video  link.  In  flight  simula¬ 
tors,  tbe  electronic  map  can  help  reduce  the  equipment  cost  of  genera¬ 
ting  terrain  displays.  Exploration  of  these  and  other  applications  has 
been  Initiated. 

In  terms  of  circuit  simplicity,  speed,  and  cost,  the  AETMS  achieves 
near-optimal  efficiency.  Although  discussion  has  focused  on  avionics 
applications,  many  other  systems  can  profitably  utilize  the  technology 
described  herein. 
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Figure  2.  Memory  Hierarchy  and  Data  Flow 
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Figure  3.  Airborne  Terrain  May  and  Display  System 
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(U)  ASSESSMENT  OF  GaAs  PASSIVATION  AND  ITS  APPLICATIONS 

Abstract 

Within  the  Air  Force,  there  is  a  need  for  a  digital  pro¬ 
cessing  capability  requiring  clock  rates  greatly  exceeding 
those  possible  with  even  the  most  advanced  silicon  technology. 
Projected  1980  to  1985  requirements  for  electronic  warfare, 
telemetry  and  digital  communications  cover  the  1  to  10  GHz 
clock  frequency  range.  Progress  is  being  made  in  fabricating 
GaAs  MESFETs  for  analog  applications  and  in  developing  GaAs 
digital  MESFET  integrated  circuits. 

One  area  of  great  importance  to  this  technology  is  in 
passivating  the  surface  of  the  GaAs  structure  in  order  to: 

(a)  mechanically  and  chemically  protect  the  device,  (b)  elec¬ 
trically  control  generation  and  recombination  of  charge  car¬ 
riers  at  the  surface,  and  (c)  improve  the  breakdown  voltage 
and  the  leakage  current  of  diodes.  Whereas  in  the  silicon 
technology,  Si02  provides  an  excellent  passivation  layer,  no 
adequate  passivation  layer  exists  presently  for  GaAs  devices. 
Consequently,  surface  passivation  of  GaAs  is  an  area  of  active 
research. 

AFAL  has  been  working  in-house  on  GaAs  passivation  using 
wet  anodized  films.  These  films  are  being  investigated  chemi¬ 
cally  with  Auger  and  ESCA  techniques.  Electrical  measurements 
are  performed  to  evaluate  capaci tance- vol tage  characteristics. 
A  transient  photo-charge  measurement  technique  was  developed 
to  evaluate  band  bending  at  the  semiconductor  surface  and  to 
determine  the  oxide  capacitance.  These  measurements  also 
allow  quantitative  determination  or  charge  injection  into  the 
insulator. 

A  major  effort  will  be  devoted  to  the  application  of  this 
passivation  technology  to  active  GaAs  devices,  especially  to 
GaAs  MESFETs. 
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I. 


INTRODUCTION: 


Today's  military  electronic  systems  use  Si  integrated 
circuits  as  the  predominate  active  component  due  to  the  ad¬ 
vanced  status  of  development  of  this  technology.  However,  the 
speed  and  the  frequency  of  operation. of  these  devices  is 
limited  (Figure  1).  There  exist  many  applications  within  Air 
Force  systems,  where  higher  speed  or  higher  frequency  of 
operation  of  electronic  circuits  would  significantly  improve 
the  performance  of  those  systems  (Figure  2).  One  example  is: 
radar  signal  processing,  where  higher  resolution  is  proportional 
to  processing  speed.  Consequently,  DOD  and  especially  AFAL  are 
developing  GaAs  technology  for  high  speed  applications. 

GaAs  is  a  semiconductor  much  like  silicon.  However,  the 
electron  mobility  of  GaAs  is  approximately  5000  and  whereas  Si 
is  only  1000  cm^/volt  sec,  or  in  other  words  the  GaAs  mobility 
is  approximately  five  times  that  of  Si  (Figure  3)1.  This 
advantage  in  electron  mobility  is  directly  reflected  in  the 
speed  of  operation  of  GaAs  devices.  There  exist  many  different 
semiconductors  with  equal  or  even  better  electron  mobilities 
(Figure  4).  However,  GaAs  technology  has  received  a  large 
amount  of  development  for  different  applications,  such  as: 
lasers,  detectors,  LED's,  and  microwave  devices.  For  these 
reasons,  much  interest  and  funding  are  expended  in  further 
developing  GaAs  technology  and  it  is,  therefore,  much  further 
ahead  of  other  compound  semiconductors. 

Due  to  the  different  material  properties  of  Si  and  GaAs, 
different  types  of  transistors  are  utilized  (Figure  5).  In  Si 
technology,  mainly  bipolar  transistors  and  insulated  gate 
field  effect  transistors  are  utilized.  Due  to  the  short  life¬ 
time  of  minority  carriers  (in  the  order  of  picoseconds),  bi¬ 
polar  transistors  cannot  be  built  in  GaAs.  Furthermore,  no 
insulated  gate  field  effect  transistors  of  operational  quality 
are  being  fabricated  due  to  limitations  in  the  insulator.  The 
main  technology  of  GaAs  circuits  is  based  on  Schottky  barriers 
or  Metal -Semi conductor-Fi el d-Effect-Transi stors  (MESFETS). 

As  indicated  in  Figure  5,  MESFETS  are  formed  utilizing 
semiconducting  GaAs  islands  on  semi-insulating  GaAs  substrates. 
The  islands  are  formed  either  by  epitaxial  growth  or  by  ion 
implantation  of  suitable  dopants  into  semi- i nsul ati ng  GaAs. 

In  operation,  the  gate  voltage  controls  the  conductivity  of 
the  gate  area.  Typically,  n-type  GaAs  islands  are  used,  i.e.: 
the  mobile  charge  carriers  are  electrons.  With  zero  volts 
applied  to  the  gate  with  respect  to  the  source,  the  transistor 
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conducts  current  between  drain  and  soun^ce.  However,  when  a 
negative  potential  of  minus  three  volts  is  applied  to  the 
gate,  the  gate  area  is  depleted  or  in  other  words,  the  elec¬ 
trons  are  repelled  from  the  gate  area.  Under  this  condition, 
no  conductance  exists  between  drain  and  source.  This  example 
illustrates  how  the  drain  to  source  conductance  can  be  con¬ 
trolled  by  the  gate  voltage.  The  gate  represents  a  Schottky 
diode  which  is  blocking  with  negative  applied  voltages;  hence, 
no  gate  current  flows  under  normal  operation  of  the  device. 

Considerable  work  is  in  progress  to  build  integrated  cir¬ 
cuits  in  GaAs,  utilizing  MESFET  technology.  At  present, 
medium  scale  integrated  circuits  are  being  fabricated,  consist¬ 
ing  of  up  to  100  gates.  A  typical  circuit  is  a  "divide-by¬ 
eight  counter"  (Figure  6),  which  operates  with  input  frequencies 
from  d.c.  to  8  GHz.^  Work  is  in  progress  to  improve  manufac¬ 
turability,  reliability,  performance  and  density  of  these 
ci rcui ts . 

One  main  problem  with  the  GaAs  technology  is  the  absence 
of  a  suitable  insulator  technology  for  surface  passivation. 

The  main  reason  for  the  fast  development  of  Si  integrated  cir¬ 
cuits  lays  on  the  existence  of  an  almost  ideal  Si  based  insul¬ 
ator  technology,  namely:  Silicon  Dioxide  (Si02).  This  insul¬ 
ator  can  be  grown  by  thermal  oxidation  of  the  5i  substrate  in 
an  oxidizing  atmosphere,  such  as  dry  oxygen  (Figure  7).  The 
Si-SiOp  interfaces  can  be  prepared  to  be  almost  perfect,  i.e.; 
they  contain  only  minute  (10^0/cn^2)  surface  state  densities 
(Figure  8).  These  insulators  can  be  advantageously  employed 
in  insulated  gate  field  effect  transistors  and  charge  coupled 
devices.  Furthermore,  this  insulator  can  be  used  to  passivate 
Si  surfaces  (Figure  9).  One  example  is  the  passivation  of 
surfaces  over  p-n  junctions  to  minimize  surface  leakage  cur¬ 
rents  and  to  maximize  breakdown  voltages.  Also,  Si02  can  be 
utilized  to  passivate  Si  surfaces  from  chemical  attack  due  to 
gases  and  vapors.  Finally,  SiOo  can  be  used  as  a  doping 
barrier  during  doping  of  special  areas  of  the  device,  such  as: 
source  and  drain  in  a  MOSFET. 

Much  work  has  been  devoted  towards  establishment  of  an 
acceptable  insulator  for  GaAs  devices.  From  these  studies  we 
realize,  that  it  is  very  unlikely  that  an  universal  insulator 
technology  will  be  found  for  GaAs  which  will  be  similar  in  all 
or  most  respects  to  Si02  on  Si  (Figure  10).  Hence  work  is 
underway  to  develop  specific  insulators  for  specific  requirements 
in  GaAs  device  technology.  One  area  of  critical  importance 
is  the  passivation  of  the  area  between  drain  and  gate  and 
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between  source  and  gate  in  GaAs  MESFETs.  The  surfaces  contain 
large  surface  state  densities  which  change  their  charge  with 
the  change  in  potential  in  the  specific  area.  This  charging 
results  in  a  reduction  of  amplification  in  these  transistors, 
especially  at  high  frequencies. 

Consequently,  AFAL  is  working  on  developing  a  passivation 
technology  for  GaAs  MESFET  application.  Due  to  the  complexity 
of  the  problem,  a  number  of  insulator  technologies  have  been 
evaluated  and  new  analytical  techniques  have  been  studied 
which  are  described  in  the  following  sections. 

II.  INSULATOR  FORMATION  TECHNIQUES: 

1 .  Deposited  Films: 

Deposited  films  of  silicon  oxide,  silicon  nitride, 
silicon  oxynitride  and  germanium  nitride  have  all  been  studied. 

In  these  experiments,  a  main  concern  is  the  compatibility  of 
the  insulating  film  with  the  GaAs  substrate.  First  of  all, 
the  adherence  of  the  film  to  the  substrate  must  be  sufficient 
to  allow  processing  of  the  structure  without  liftoff  of  the 
insulating  film.  Furthermore,  the  thermal  expansion  coefficients 
must  be  sufficiently  matched;  problems  in  this  area  cause 
breakage  of  the  GaAs  substrate  since  GaAs  is  more  brittle  than 
Si.  In  this  respect,  it  is  desirable  to  deposit  these  films 
at  not  too  high  temperatures  (say  ZSC^C  instead  of  750°C)  to 
minimize  damaging  effects  due  to  mismatch  in  the  expansion  co¬ 
efficients.  Advanced  deposition  techniques,  such  as:  plasma- 
enhanced  deposition  or  photo-enhanced  deposition  techniques 
are  being  used.  Of  main  interest  is  the  establishment  of  an 
insulator-semiconductor  interface  with  sufficiently  low  sur¬ 
face  state  density. 

2 .  Thermal  Oxidation: 

Another  approach  being  actively  pursued  is  the 
formation  of  a  thermal  oxide  but  not  necessarily  gallium  or 
arsenic  oxide.  The  oxidation  of  GaAs  requires  temperatures  of 
approximately  500°C.  At  these  temperatures,  As^O,  is  vola¬ 
tile,  leaving  Ga^Coj  which  appears  as  a  polycrystalline, 
conductive  film  rather  than  an  amorphous,  insulating  film. 

The  reason  for  this  is  due  to  the  fact  that  Ga  is  a  conditional 
rather  than  a  primary  glass  former  and  requires  another  ele¬ 
ment  to  form  gl assy  materi al .  Experiments  on  thermal  oxida¬ 
tion  of  GaAs  Pj^_^indi cate  that  the  thermally  formed  oxides  are 
glassy  and  i  fiiSu  fati  ng .  One  approach  therefore,  is  to  form  a 
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GaiAs  P-.  surface  layer  on  GaAs  by  1on  1  pi  antatlon  ot  p  and 
subsequently  to  thermally  oxidize  this  structure.  The  advan¬ 
tage  of  this  approach  is  that  the  semiconductor-insulator 
interface  is  formed  within  the  GaAs  substrate  and  consequently, 
this  interface  was  at  no  time  exposed  and/or  contaminated  by 
the  atmosphere.  In  contrast,  one  of  the  main  concerns  with 
deposited  films  is  that  the  semiconductor  surfaces  become  the 
interfaces  and  consequently  are  contaminated  to  some  extent. 

A  further  approach  for  formation  of  insulating  films 
utilizes  the  anodic  process;  at  present,  most  of  our  insulating 
films  are  prepared  by  wet  anodization.'^ 

3 .  Aqueous  Anodic  Oxidation: 

When  the  semiconductor  is  inserted  into  an  aqueous 
solution  with  its  high  dielectric  constant,  the  surface  atoms 
of  the  semiconducting  material  are  dissolved  into  the  liquid 
as  their  ionic  bonds  are  weakened.  The  atoms  enter  the  liquid 
as^tpi^l  ti  pi  e  positively  ionized  species,  such  as  Ga"^"^  or 
As  ,  leaving  on  the  semiconductor  surface  an  immobile  nega¬ 
tive  charge  which  has  to  be  neutralized  by  mobile  holes  from 
the  bulk  of  the  semiconductor.  After  a  second  electrode  is 
placed  as  cathode  electrode  into  the  solution,  the  semiconduc¬ 
tor  can  be  employed  as  anode  with  a  positive  potential  which 
aids  the  hole  flow  to  the  semiconductor  surface  so  that  more 
ions  can  be  released.  These  ions  travel  towards  the  cathode. 

If  sufficient  ions  are  released  they  form  a  dense 
ionic  cloud  in  front  of  the  anode.  This  is  particularly  easy 
if  either  the  solubility  or  the  diffusion  coefficient  of  these 
ions  is  reduced  by  the  introduction  of  a  suitable  liquid  into 
the  aqueous  solution.  It  was  found  that  some  liquids  are  par¬ 
ticularly  effective,  such  as  glycol. 

The  ions  in  this  cloud  will  capture  OH"  ions  of  the 
aqueous  solution,  and  with  the  simultaneous  formation  of  HgO, 
are  oxidized.  It  is  likely  that  most  of  them  will  form  weTl- 
known  oxides  such  as  GajOo  or  As^Ogj  which  will  exceed  a 
density  saturation  lever  Tn  this^cToud  so  that  they  begin  to 
precipitate  against  the  semiconductor.  An  oxide  layer  begins 
to  get  deposited  thus. 

Further  growth  is  now  only  possible  if  ionic  trans¬ 
port  takes  palce  across  the  oxide.  Each  of  the  ionic  species 
has  a  threshold  field  for  drift  and  a  specific  value  of  mobility. 
With  normal  current  densities  applied  across  the  semiconductor 
surface  which  is  immersed  into  the  electrolyte,  the  fields 
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across  the  oxide  and  its  interfaces  are  sufficiently  high 
for  the  motion  of  all  ions,  namely  OH'  and  the  positively 
ionized  components  of  the  semiconductor.  However  as  the 
mobility  of  OH'  is  much  smaller  than  that  of  the  multiply 
ionized  components,  the  primary  growth  occurs  at  or  near  the 
electrolyte-oxide  interface.  Unfortunately,  the  current 
density  cannot  be  reduced  beyond  a  minimum  value  which  is 
different  for  each  material.  This  value  has  sometimes  been 
called  the  dissolution  current  as  it  seems  to  agree  with  the 
equivalent  cation  current  due  to  the  dissolution  of  the  oxide 
in  the  electrolyte  when  no  external  current  is  applied.  It 
can  be  understood  as  the  flow  of  those  cations  which,  after 
dissolution  from  the  oxide,  successfully  leave  the  anode 
region.  Only  when  this  minimum  current  density  is  exceeded, 
can  growth  take  place.  This  means  that  with  the  electrolytes 
known  so  far,  native  oxides  on  most  III-V  Compound  Semicon¬ 
ductors  such  as  GaAs,  can  only  be  produced  after  exceeding  the 
drift  threshold  for  all  the  ionic  species  involved. 

However,  A1  and  some  other  metals  have  a  Very  small 
dissolution  current.  After  depositing  a  layer  of  such  a  metal 
by  evaporation  on  top  of  the  semiconductor,  this  metal  will 
first  be  oxidized  in  the  electrolyte.  After  it  has  been  com¬ 
pletely  transformed,  the  oxidation  of  the  semiconductor  be¬ 
gins.  For  a  high  current  density,  the  semiconductor  components 
drift  as  cations  across  the  AKOo  and  native  oxide  of  the 
semiconductor  appears  to  top  of  the  AKO^.  On  the  other  hand, 
for  a  low  current  density  the  semi conducfor  components  cannot 
become  mobile,  and  it  is  only  OH”  which  crosses  the  Al^O,. 
Correspondingly,  native  oxide  of  the  semiconductor  appears 
only  between  AlgO,  and  the  semiconductor.  It  is  therefore 
possible  to  produce  various  types  of  layered  structures.  It 
is  useful  to  consider  this  now  systematically  for  various 
types  of  device  applications. 

The  oxide  growth  can  be  undertaken  either  with  a 
constant-current  source  or  with  a  constant-voltage  supply  in 
series  with  a  sufficiently  high  series  resistance  in  order  to 
limit  the  initial  current  density.  A  voltage  supply  without 
such  a  limiting  resistance  is  not  advisable  as  the  initial 
current  density  can  easily  become  so  high  that  hot  spots  and 
other  detrimental  effects  can  occur  during  the  initial  stage 
of  growth  so  that  a  deterioration  of  the  oxide  quality  appears. 
A  constant-current  case  would  require  a  voltage  measurement 
between  cathode  and  anode  in  order  to  monitor  the  growth  pro¬ 
cedure.  A  voltage  and  series-resistance  agreement  can  enable 
one  to  monitor  the  growth  by  measuring  the  current  flowing 


across  the  electrolytic  cell.  For  the  relatively  high  final 
growth  voltages  involved  across  most  of  the  oxides  of  device 
interest,  it  is  not  required  to  undertake  the  more  refined 
electrochemical  potential  measurements  as  these  would  then  be 
normally  a  negligible  fraction  of  the  applied  voltage. 

It  is  possible  to  model  the  growth  behavior  rela¬ 
tively  satisfactorily  by  a  simple  treatment,  provided  that 
precautions  are  taken  for  a  uniform  current  density  across  the 
semiconductor  surface  which  is  exposed  to  the  electrolyte.  By 
way  of  example,  this  is  to  be  outlined  here  for  the  constant- 
voltage  plus  series-resistance  arrangement.  The  treatment  of 
the  other  case  of  growth  is  very  similar. 

The  general  current  functions  measured  are  shown  in 
Figure  11.^  A  typical  experimental  set  up  is  given  in  Figure 
12,  where  it  is  also  indicated  that  a  light  source  is  needed 
in  order  to  have  sufficient  holes  available  for  the  neutrali¬ 
zation  of  the  negative  surface  charge  created  by  the  removal 
of  positive  cations.  This  illumination  is  particularly  important 
for  n-type  material  where  practically  no  holes  would  be  present 
otherwise.  Failing  to  provide  a  sufficient  supply  of  holes 
means  that  the  space  charge  layer  set  up  by  the  negative 
surface  charge  caused  by  the  departing  cations  absorbs  a 
substantial  part,  if  not  all,  of  the  externally  applied 
voltage  so  that  the  oxide  thickness  will  be  much  smaller  than 
that  predicted  by  the  model. 

Recent  experiments  were  performed  also  utilizing 
non-aqueous  anodic  growth  (rather  than  the  above  described 
aqueous  anodi c  growth) ,  utilizing  an  ammonium  pentaborate 
solution  in  ethylene  glycol.  Superior  results  were  obtained 
when  A1  films  on  GaAs  were  anodized  in  this  solution. 

4 .  Plasma  Anodization: 

We  are  now  in  the  process  of  establishing  a  plasma 
anodization  growth  facility.  The  plasma  anodization  will 
result  in  less  contaminated  films  than  the  wet  anodized  films 
since  gases  are  obtainable  with  less  contamination  content 
than  liquids.  The  attached  Figure  13  indicates  the  setup  for 
plasma  anodic  growth. 6  An  R.F.  discharge  creates  an  oxygen 
plasma  and  oxygen  ions  or  excited  neutral  oxygen  atoms  reach 
the  surface  of  the  sample  to  be  anodized.  A  voltage  applied 
to  the  sample  creates  an  electric  field  in  the  growing  insul¬ 
ator  of  approximately  10®  V/cm  which  assists  in  the  migration 
of  the  ions,  namely  0,  Ga,  As,  and  A1 .  Superior  films  have 
been  reported  when  A1  films  on  GaAs  were  formed. 


III.  ANALYTICAL  TECHNIQUES: 


The  insulating  films  on  GaAs  and  their  interfaces  are 
evaluated  analytically  for  their  mechanical  stability,  for 
their  chemical  characteristics  and  composition  and  for  their 
electrical  behavior. 

The  insulating  films  have  to  have  adequate  adherence  to 
the  substrate.  The  films  have  to  appear  uniform  on  inspection 
with  a  microscope,  and  liftoff,  peeling  of  the  film  or  rupture 
of  the  substrate  is  not  acceptable,  since  the  films  have  to 
stand  further  . process! ng  such  as  heat  treatment  and  immersion 
into  liquids,  care  must  be  taken  that  these  processes  don't 
damage  the  films. 

In  recent  years,  instrumental  techniques  for  chemical 
analysis  of  thin  layers  have  been  greatly  improved.  A  survey 
of  these  techniques  was  recently  presented  by  G.  R.  Larrabee.^ 
We  have  utilized  optical  characterization.  X-ray  Photoelectron 
Spectroscopy  (or  ESCA)>  Secondary  Ion  Mass  Spectroscopy 
(SIMS)  and  Auger  Electron  Spectroscopy  (AES).  To  detail  these 
techniques  and  describe  all  of  our  results  would  be  beyond  the 
scope  of  this  presentation.  However,  to  give  an  example  of 
the  work  involved,  we  describe  here  the  AES  experiment  per¬ 
formed  in  an  ultra-high  vacuum  system  (Fig  14).  The  sample  is 
bombarded  with  electrons  of  several  thousand  electron  volts 
(keV)  energy.  These  high  energy  electrons  are  able  to  ionize 
atoms  in  the  sample  by  knocking  out  electrons  from  inner 
shells.  Now  electrons  with  higher  energy  can  relax  into  the 
empty  state,  thereby  releasing  a  characteristic  energy.  This 
energy  is  transferred  to  another  electron  which  subsequently 
is  emitted  into  vacuum  where  it  is  detected  and  its  energy 
measured  and  related  to  specific  atoms.  The  attached  Figure 
15  indicates  the  results  obtained  showing  the  derivative  dN/dE 
signal  vs.  energy  E.  Two  different  traces  are  shown  indica¬ 
ting  different  surface  conditions.  Simultaneous  to  the  elec¬ 
tron  bombardment,  the  samples  are  bombarded  with  argon  ions 
with  energies  of  several  keV.  These  argon  ions  sputter  off 
the  surface  of  the  sample  thereby  profiling  the  film  composi¬ 
tion.  Figure  16  presents  such  a  profile  of  an  anodized  film 
of  A1  on  GaAs.  On  the  abscissa  the  sputter  time  is  presented 
which  relates  to  the  distance  from  the  film  surface  and  on  the 
ordinate,  the  AES  signal  strength  for  the  important  species 
are  presented  namely  0,  A1 ,  Ga  and  As.  These  data  let  us 
understand  the  chemical  composition  of  the  film  at  different 
locations  and  especially  allow  us  to  evaluate  the  chemical 
composition  in  the  insulator-semiconductor  transition  region. 
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The  most  important  studies  evaluate  the  electronic 
behavior  of  these  films  and  their  interfaces  since  adequate 
electronic  characteristics  are  required  for  proper  device 
performance.  The  most  common  technique  is  to  measure  capa¬ 
citance  as  a  function  of  voltage  (C-V)  for  metal -insulator- 
semi  conductor  structures.^  In  general  the  sample  is  first 
biased  to  attract  semiconductor  charges  to  the  interface  or  to 
put  the  interface  into  accumulation  (positive  voltage  for  n- 
type  semiconductor).  The  capacitance  measured  then  is  that  of  the 
insulator.  Now  the  voltage  is  altered  to  a  reverse  position 
(negative  voltage  for  n-type  semiconductor),  causing  depletion 
of  carriers  under  the  metal  electrode.  Hence  the  insulating 
structure  now  consists  of  the  depleted  semiconductor  and  the 
insulating  oxide  film,  thereby  decreasing  the  capacitance. 

Figure  17  demonstrates  a  typical  C-V  characteristic  for 
anodized  GaAs  structures.  Difficulties  have  been  encountered 
in  interpreting  these  capacitance  characteristics,  since 
charge  apparently  is  injected  into  the  insulator  and  this 
charge  influences  the  capacitance  values  measured.  In  order 
to  better  evaluate  these  measurements,  we  incorporated  a 
photopulse  technique.  In  a  pulse  setup,  we  first  put  the 
metal-insulator-semiconductor  structure  into  depletion  and 
measure  within  10  milliseconds  the  capacitance  of  the  struc¬ 
ture;  we  then  shine  a  strong  photopulse  of  approximately  4 
milliseconds  duration  onto  the  sample,  which  causes  the  deple¬ 
tion  region  to  collapse  and  the  generation  of  a  voltage  pulse 
of  the  size  of  the  surface  potential.  This  method  then  allows 
direct  evaluation  of  the  capacitance  as  a  function  of  surface 
potential,  thereby  eliminating  the  influence  of  the  unquan¬ 
tified  charging  of  the  insulator.  These  data  can  then  be 
analyzed  theoretically  to  obtain  information  on  doping  concen¬ 
tration.  These  transi ent-photo-measurement  results  can  then 
be  applied  to  the  C-V  results  to  evaluate  charging  of  the 
films.  We  also  measure  the  charging  and  discharging  of  the 
sample  under  the  above  conditions  to  obtain  direct  information 
on  the  capacitance  of  these  structures. 

IV.  DEVICE  APPLICATIONS: 

A  large  amount  of  our  future  work  will  be  concerned  with 
applying  this  passivation  technology  to  GaAs  MESFETs  and  GaAs 
MESFET  integrated  circuits.  Our. first  concern  will  be  to 
apply  this  technology  to  the  areas  between  gate  and  drain  and 
between  gate  and  source  of  MESFETs.  Further  use  will  be 
■  concerned  with  dielectric  isolation  of  these  MESFETs  to  simplify 
processing  and  to  improve  density  and  performance  of  these 
circuits.  We  will  keep  an  open  mind  to  consider  at  intervals 
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the  applicability  o£  these  passivation  layers  to  depletion 
type  insulated  gate  FETs  as  well  as  to  GaAs  Charge  Coupled 
Devices  (CCDs) . 

Due  to  the  very  high  speed  of  GaAs  MESFET  circuits 
primary  applications  will  be  in  very  high  speed  logic 
(i.e.:  10-18  GB/S)  for  radar  signal  processing,  ECM, 
satellite  communications  and  computer  applications.  A  major 
impact  will  be  from  the  fabrication  of  very  high  speed  A/D 
converters.  Since  GaAs  integrated  circuits  will  be  premium 
cost  circuits,  their  use  will  be  restricted  to  the  highest 
speed  parts  of  the  system  or  to  special  applications  where 
silicon  circuits  will  not  function. 

If  the  development  of  this  insulator  technology  progresses 
far  enough,  MOSFET  type  devices  may  become  feasible  in  GaAs 
and  these  could  be  produced  at  lower  cost  due  to  their  higher 
density  and  lower  power  dissipation.  It  is  unliTcely,  however, 
that  these  GaAs  MOSFET  integrated  circuits  would  ever  be 
cost  competitive  with  silicon  circuits  at  speeds  where 
silicon  can  be  used. 

Transferred  electron  devices  are  essentially  ultra- 
high  speed  switching  devices  which  operate  at  picosecond 
speeds  with  power  dissipations  of  about  60  mW/gate  and  very 
high  field  strengths  to  initiate  a  "high  field  domain" 
which  propagates  down  the  device  at  about  10^  cm/sec. 

Gate  areas  are  typically  about  20  mil^  whereasTTL  gates  are 
typically  about  110  mil'^,  S^hottky  TTL  gates  about  126  mil'^, 
and  ECL  gates  about  120  mil^.  The  transferred  electron  effect 
occurs  in  several  III-V  semiconductor  materials  (i.e.:  GaAs 
and  InP)  due  to  a  bulk  negative  differential  mobility  caused 
by  the  transfer  of  electrons  from  a  high  mobility  band  to 
a  low  mobility  band.  Since  TEDs  are  high  field  devices, 
achievement  of  high  fields  with  low  leakage  demands  high 
quality  insulators  for  passivation  of  the  surface  in  planar 
processes.  Likewise,  the  achievement  of  low  noise  at  high 
fields  is  critically  dependent  on  low  leakage  insulators  and 
passivation. 

One  o£  the  major  applications  of  GaAs  TEDs  will  be  in 
A/D  converters.  The  performance  of  analog  to  digital  converters 
for  ELINT  applications  is  limited  by  the  sampling  aperture 
width.  GaAs  TEDs  will  permit  a  5  picosecond  sampling  aperture 
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with  TED  comparators  which  would  allow  frequency  components 
up  to  500  MHz  with  6  bit  resolution.  Competitive  tunnel  diode 
sample  and  hold  circuits  are  limited  to  about  16  pS  or 
greater  aperture.  TEDs  provide  a  220%  improvement  in  band¬ 
width. A  noise  threshold  of  50-100  mV  requires  each  level  of 
an  A/D  converter  to  be  200  mV  or  more,  therefore,  for  63 
levels,  this  means  at  least  12.6  V  input  is  required  into 
the  digital  converter.  This  in  turn  requires  high  voltage 
amplifiers  with  wide  bandwidth  (which  are  not  currently 
available) .  Suitable  passivation  should  be  able  to  cut  down 
this  noise  level  and  make  these  devices  much  more  practical. 

In  general,  the  most  fruitful  applications  for  GaAs  TEDs 
are  probably  in  the  following  areas:  A/D  and  D/A  converters, 
frequency  dividers,  latches,  special  purpose  gates,  shift 
register  memories,  counters,  and  decoders. 

GaAs  CCDs  show  great  potential  for  use  in  high  speed 
image  sensors  in  the  IR  and  visible  regions  of  the  spectrum. 
Since  good  detector  devices  can  be  built  in  GaAs  -  the  problem 
of  high  speed  image  scanning  could  be  solved  with  a  minimum 
of  hardware,  interconnects,  and  with  high  reliability  at 

low  cost -  if  an  adequate  CCD  device  with  high  transfer 

efficiency  and  low  leakage  can  be  fabricated.  The  dielectric 
films  for  isolation  and  passivation  of  this  device  are  a  key 
technology  for  the  success  of  this  new  device.  Planar, 

Schottky  gate,  buried  channel  CCDs  require  very  narrow  gaps 
between  the  metal  Schottky  gates  resulting  in  a  configuration 
with  relatively  large  areas  of  metal  separated  by  narrow 
gaps  and  substantially  long  lengths  of  the  gap  per  transfer 
stage.  This  configuration  contrasts  with  MESFETs  where 
narrow  stripes  of  metal  are  surrounded  by  large  clear  areas. 
Two  level  metallizations  with  a  large  number  of  cross  overs 
per  device  require  a  deposited  insulator  with  very  low  pin¬ 
hole  density  and  excellent  step  coverage  as  well  as  good 
adhesion  and  dielectric  strength.  Currently. a  plasma-enhanced, 
deposited  glassy  dielectric  of  proprietary  composition  is 
being  used  to  fabricate  early  experimental  GaAs  CCDs  (Figure 
18). 

In  summary,  the  passivation  of  GaAs  surfaces  with  a 
suitable  high  quality  insulator  is  necessary  to  achieve  low 
surface  state  density,  low  leakage,  and  low  noise  in  a 
number  of  advanced  new  devices  (such  as:MESFETs,  MOSFETs , 

TEDs,  and  CCDs).  Air  Force  applications  of  these  new  devices 
will  be  primarily  in  premium  cost,  high  speed  circuits  where 
performance  at  above  10  GB/S  is  of  paramount  importance. 

Many  such  applications  exist  within  the  Air  Force  in  the.', 
areas  of  radar,  ECM,  satellite  communications,  high  speed 
computers  and  image  sensors.  Application  of  this  advanced 
technology  can  reduce  the  cost  of  critical  parts  of  systems 
where  a  large  premium  is  paid  for  high  performance. 
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SPEED  OF  CURRENT  Si  INTEGRATED  CIRCUITS 


POWER  DISSIPATION  (MW/GATE) 


AIR  FORCE  APPLICATIONS  FOR 
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Figure 
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MOBILITY  OF  VARIOUS  SEMICONDUCTORS 
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GaAs  DIVIDE  BY  8  COUNTER 
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PLASMA  ANODIZATION  SYSTEM 
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AUGER  ELECTRON  SPECTROMETER 
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PROFILE  OF  ANODIZED  A1  ON  GaAs 


CAPACITANCE  -  VOLTAGE  CHARACTERISTICS  OF  MOS  GaAs 
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GaAs  SCHOTTKY  BARRIER  CCD 


TRANSFERRED  ELECTRON  DEVICE 


Figure  19 
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Identification  of  Impurity  Complexes  in  Gall i um  Arsenide 
Device  Material  by  High  Resolution  Magneto-PhOtol umi nescence 

Abstract 

Advanced  avionics  systems  under  development  or  proposed 
for  future  applications  rely  heavily  on  the  capabilities  of 
gallium  arsenide  (GaAs)  as-  an  electromagnetic  material.  Many 
of  the  characterization  techniques  applicable  to  silicon  have 
required  considerable  refinement  for  GaAs  and  other  more  com¬ 
plicated  compound  semiconductors.  One  such  technique  is 
laser-pumped  photoluminescence  with  and  without  applied 
external  magnetic  and  stress  fields.  This  measurement  in¬ 
volves  exciting  electrons  in  the  material  using  low  power 
laser  radiation  and  observing  the  light  emitted  by  the  elec¬ 
trons  as  they  relax.  This  emission  can  be  used  to  characterize 
impurities  and  defects  since  both  create  distinguishable 
electronic  states.  Many  impurity  states  in  GaAs  have  been 
investigated  in  depth  and  are  well  understood. 

This  study  concerns  five  less  understood  transitions 
which  we  have  investigated.  These  systems  are  not  identifi¬ 
able  with  any  simple  impurity  transitions  having  similar 
binding  energies.  When  perturbing  fields  are  applied,  dif¬ 
ferences  in  the  behavior  of  these  transitions  and  simple 
impurity  transitions  are  observed.  Based  on  these  effects, 
models  have  been  developed  to  explain  the  sources  of  these 
transitions.  The  model  currently  proposed  involves  a  compli¬ 
cated  molecule-like  arrangement  of  an  anti-site  defect  and 
nearby  impurities.  This  measurement  is  the  first  experimental 
verification  of  the  existence  of  such  defects. 

The  development  of  high  frequency  or  high  power  devices 
could  be  limited  by  these  i nherent  defects  which  cannot  be 
eliminated  by  further  purification  processes.  Our  research 
effort  is  motivated  by  the  need  to  understand  and  eliminate 
or  control  such  defects  in  the  demanding  environment  of 
advanced  Air  Force  avionics  devices. 
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Introduction 

The  feasibility  and  economic  attractiveness  of  many 
future  Air  Force  Avionics  and  ground-based  electronics  sys¬ 
tems  require  the  development  of  new,  smaller,  lighter,  more 
inexpensive,  and  longer  lived  electronic  devices.  Other  sys¬ 
tems  ,  less  critically  dependent  on  the  development  of  new  de¬ 
vices,  will  derive  tremendous  cost  and  performance  benefits 
from  such  developments.  Phased  array  radar  systems,  for 
example,  offer  vast  improvements  in  performance  over  conven¬ 
tional  radar  systems,  but  their  implementation  is  generally 
limited  to  frequencies  below  about  8  GHz  primarily  because  of 
a  lack  of  suitable  microwave  power  devices.  Developments 
leading  to  improved  power  and  efficiency  of  devices  for  use 
in  phased  array  radar  systems  would  also  have  a  revolutionary 
impact  on  other  systems,  such  as  satellite  communication 
transmitters,  for  example.  Similar  developments  would  also 
address  the  far-reaching  requirements  for  electromagnetic 
warfare  systems,  whose  designs  encompass  digital  information 
processing  speeds  currently  impossible  to  achieve  and  even 
beyond  the  projected  optimum  capabilities  of  silicon  tech¬ 
nology. 

In  each  of  these  specific  systems,  it  has  been  the  po¬ 
tential  performance  of  gallium  arsenide  devices  which  has 
permitted  the  prediction  of  solutions  for  critical  design 
problems.  Laboratory  feasibility  experiments  bn  the  perfor¬ 
mance  of  GaAs  field  effect  transistors  have  indicated  that 
this  device  is  one  likely  solution  to  technological  roadblocks 
limiting  incorporation  into  our  inventory  of  many  of  the 
advanced  electronic  systems  discussed.  Such  experimentation 
has  also  demonstrated  that  the  ultimate  performance  of  FET's, 
and  of  other  new  devices  being  developed,  depends  directly  on 
the  success  with  which  materials  can  be  produced  in  a  care¬ 
fully  controlled  and  reproducible  manner. 

Materials  for  Device  Technology 

The  establishment  of  a  sound  technology  base  allowing 
reliable  and  reproducible  material  synthesis  requires  a  basic 
understanding  of  the  material  physics,  including  the  effects 
of  foreign  impurities,  host  defects  and  complexes  involving 
both  impurities  and  defects.  Native  defects,  or  crystalline 
lattice  imperfections,  can  interact  with  intentionally 
introduced  dopants  to  form  neutral  centers.  It  is  also 
possible  for  unintentional  residual  acceptors  to  compensate 
intentionally  added  donors.  In  either  case,  the  resulting 
electrical  performance  of  the  material  will  differ  from  the 
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desired  design  parameters  and  will  adversely  affect  device 
performance. 

A  major  thrust  of  our  research  effort  is  the  charac¬ 
terization  of  such  residual  compensating  centers,  either  de¬ 
fects  or  impurities.  Great  care  both  in  sample  growth  and 
preparation  and  in  the  application  of  sophisticated  measure¬ 
ment  techniques  is  required  to  complete  this  thrust. 


The  ef f ecti ve-mass- 1 i ke  donors  in  Ill-V  binary  and  in 
several  of  the  III-V  ternary  systems  are  very  weakly  bound. 

The  chemical  shifts  and  central  cell  corrections  which  allow 
various  impurities  and  defects  to  be  differentiated  in  terms 
of  their  characteristic  energy  separations  are  small.  This 
requires  low  residual  concentrations  to  prevent  concentration 
broadening  and  merging  of  the  impurity  levels.  Sophisticated 
crystal  growth  procedures  are  required  to  control  the  residual 


impurities  and  defects  at  an  acceptably  low  level  to  permit 
effective  controlled  doping.  Parameters,  such  as  the  purity 
of  growth  materials,  temperature  gradients  within  the  growth 
regions,  growth  rates,  dopant  diffusion  coefficients,  and 
perhaps  many  not  yet  identified,  affect  the  success  of  such 
crystal  growth  efforts.  Once  the  residual  impurities  and 
defects  have  been  controlled  at  an  acceptably  low  level, 
dopants  can  be  added  to  the  growing  system  intenti oni al ly  to 
obtain  uncompensated  and  uniformly  doped  material.  Con¬ 


trolled  doping  experiments  with  known  donors  must  also  be  in 
the  low  concentration  range  (Nq  <5  x  lOl^cm"^)  to  permit 
identification  of  specific  donors;  this  represents  a  maximum 
of  one  dopant  atom  for  each  10^0  host  atoms. 


Our  characterization  measurements  strive  to  gain  an 
understanding  of  the  detailed  nature  of  defects,  impurities, 
and  their  interactions,  so  that  this  information  can  be  used 
by  growers  to  determine  the  optimized  conditions  for  produc¬ 
tion  of  material  with  desired  characteristics. 


Characteri zati on  Techniques 


The  electrical  properties  of  semiconductors  have  a  long 
history  of  extensive  investigation.  Many  of  the  basic  prop¬ 
erties  of  these  device  materials  were  determined  from  electri¬ 
cal  measurements.  As  the  investigations  of  semiconductors 
were  extended  to  some  of  the  larger  band  gap  materials  elec¬ 
trical  measurements  were  not  as  readily  applicable.  This, 
coupled  with  the  understanding  of  excitons  and  their  contri¬ 
bution  to  the  elucidation  of  material  properties  in  the 
1960's,  led  to  a  wide  application  of  optical  studies  to 
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semiconductor  materials.  It  was  found  that  these  materials 
reflect,  absorb,  disperse,  scatter  and  radiate  light  and  in 
general,  interact  strongly  with  the  electromagnetic  radiation 
field.  Because  of  this  strong  interaction  many  of  the  funda¬ 
mental  properties  of  these  materials,  such  as  their  energy 
band  gaps,  activation  energies  of  defects  and  foreign  impuri¬ 
ties,  effective  mass  parameters,  refractive  indicies,  dielec¬ 
tric  functions,  exciton  binding  energies  and  lattice  vibra¬ 
tion  frequencies,  can  be  determined  from  optical  experiments. 
These  experiments  cover  the  electromagnetic  spectrum,  ranging 
from  the  vacuum  ultraviolet  to  the  far  infrared.  Over  the 
past  two  decades,  optical  spectroscopy  has  been  increasingly 
etnployed  for  the  study  and  measurement  of  semiconductor 
properties  and  has  ultimately  become  a  very  powerful  experi¬ 
mental  tool. 

Intense  photoluminescence  is  observed  in  many  semiconduc¬ 
tors  at  low  temperatures.  When  analyzed  spectrally,  this 
photoluminescence  provides  an  extensive  source  of  experimental 
data  which  contributes  to  the  ultimate  identification  of  the 
electronic  states  of  impurities  and  defects  in  these  semi¬ 
conductors.  Many  sharp  lines  appear  in  such  spectra,  particu¬ 
larly  from  bound  excitons,  which  provide  a  "finger  print"  of 
the  impurities  and  defects  which  are  present  in  the  semicon¬ 
ductor  lattice. 

The  technique  of  high  resolution  optical  absorption,  re¬ 
flection  and  photoluminescence  spectroscopy  has  been  exten¬ 
sively  used  to  analyse  the  intrinsic  energy  band  parameters 
of  semiconductors,  as  well  as  their  impurity  and  defect 
states.  Intrinsic  or  free  exciton  formation  is  observed  in 
most  well  formed  crystal  structures  when  optically  excited 
with  the  proper  energy  and  at  cryogenic  temperatures.  Exam¬ 
ination  of  the  behavior  of  free  excitons  in  such  crystals  has 
led  to  a  great  deal  of  success  in  understanding  the  intrinsic 
band  structure  of  semiconductors.  Bound  excitons  have  been 
applied  equally  successfully  in  probing  the  impurity  and 
defect  structure  of  many  of  these  same  materials. 

With  the  advent  of  the  laser  as  an  exciting  source, 
optical  studies  have  become  even  more  rewarding,  particularly 
when  non-linear  effects  are  being  investigated.  They  are 
also  invaluable  as  exciting  sources  fpr  high  resolution 
photoluminescence  studies. 

Our  experimental  apparatus  for  high  resolution  spectroscopy 
is  shown  in  Fig  1.  In  these  experiments  the  sample  is  mounted 
in  a  liquid  helium  Dewar  attached  to  a  vacuum  pumping  system 
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permitting  temperatures  near  1°K  to  be  realized.  The  Dewar 
may  also  be  placed  between  the  poles  of  a  magnet  permitting 
magneto-optical  studies  to  be  made.  The  optical  source  is 
tailored  to  the  specific  materials  being  investigated.  For 
our  photoluminescence  measurements,  we  employ  a  200  mW  krypton 
ion  laser  with  principal  laser  radiation  at  647.1  nm.  The 
transmitted,  reflected,  or  emitted  light  may  then  be  analyzed 
with  a  high  resolution  spectrograph.  The  heart  of  the  spec¬ 
trograph  is  the  reflection  grating  selected  to  give  optimum 
data  for  the  particular  material  under  investigation.  The 
specific  design  parameters  for  our  spectrograph  have  been 
optimized  for  investigating  GaAs.  The  spectrograph  has  a 
four  meter  focal  length  and  employs  a  4  inch  square  grating 
with  2160  grooves/mm  blazed  at  32.7°.  This  design  gives  a 
dispersion  of  approximately  0.54A/mm  in  first  order  in  the 
free  exciton  spectral  region  of  GaAs.  Theoretical  spectral 
energy  resolution  of  0.3  cm“l  is  achievable  with  this  instru¬ 
ment.  It  may  be  employed  with  either  photographic  recording 
or  electronic  recording  of  the  optical  data. 

Photoluminescence  Spectra 

The  photoluminescence  process  in  the  intrinsic  region  is 
difficult  to  quantitatively  analyze.  In  the  case  of  absorp¬ 
tion  and  reflection  experiments  broad  band  illumination  is 
used  generally  at  relatively  low  intensities.  This  maintains 
the  crystal  in  thermal  equilibrium.  The  complications  in  the 
theory  result  from  uncertainty  in  the  proper  model  for  the 
dielectric  function  and  the  proper  surface  boundary  conditions 
These  same  problems  are  present  in  the  photol umi nescent  pro¬ 
cess  with  a  number  of  additional  problems.  Exciting  inten¬ 
sities  are  usually  greater  in  the  case  of  photoluminescence 
studies  with  laser  excitation  most  generally  used.  As  a 
result  the  crystal  may  not  be  in  thermal  equilibrium. 

Surface  recombination  as  well  as  bulk  and  surface  scat¬ 
tering  have  great  influence  on  the  photol umi nescent  results 
and  for  the  most  part  are  not  quantitatively  known.  The 
interpretation  of  intrinsic  photol umi nescent  spectra  is  in 
need  of  quantitative  theory  to  yield  as  much  reliable  infor¬ 
mation  as  is  obtained  from  absorption  and  reflection  spectra. 

These  difficulties  notwithstanding,  considerable  progress 
has  been,  and  is  being,  made  in  the  identification  of  lumines¬ 
cence  centers  in  gallium  arsenide.  The  work  described  in 
this  paper  represents  one  part  of  our  total  effort  to  develop 
photoluminescence  measurement  techniques  for  use  in  GaAs 
characterization. 
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During  the  past  four  years,  over  800  samples  of  GaAs 
have  been  characterized  in  our  facility.  We  have  focussed 
our  attention  on  the  i dentifi cation  of ,  shal  1  ow,  or  loosely 
bound,  impurities  or  defects  with  energy  levels  lying  close 
to  either  the  valence  or  conductions  bands  i n  the  materi a1 . 

These  levels  can  be  detected  when  electrons  experience  re¬ 
laxation  transitions  involving  the  levels  and  the  opposite 
energy  band.  For  example  when  ah  electron,  previously  pumped 
by  the  laser  radiation  from  an  acceptor  impurity  level  up  to 
the  conduction  band,  relaxes  back  down  to  the  impurity  level 
it  emits  a  photon  with  an  energy  equal  to  the  difference 
between  its  initial  and  final  states.  A  similar  emission  can 
be  observed  when  an  electron  relaxes  from  a  donor  level  down 
to  the  valence  band.  The  typical  emission  seen  when  this 
transition  takes  place  is  labeled  in  Figure  2.  The  line  is 
typically  broad,  which  reflects  the  energy  uncertainty  assoc¬ 
iated  with  an  electron  in  the  valence  band  of  a  semiconductor. 

A  second,  and  much  sharper,  set  of  lines  results  from  transitions 
involving  excitons  bound  to  shallow  impurities.  In  this 
context,  an  exciton  is  a  hydrogen-like  molecule  formed  by  the 
coulombic  interaction  between  a  free  electron  and  a  "hole". 

When  this  system  is  formed  in  the  presence  of  an  impurity, 
the  allowed  electron  energy  levels  become  very  sharp  with 
widths  typically  on  the  order  of  0.1cm'"'.  The  bulk  of  our 
research  effort  has  involved  the  measurement  of  these  very 
sharp-line  transitions.  Typical  donor  and  acceptor  bound 
exciton  lines  are  also  indicated  in  the  data  of  Figure  2. 

Our  work  in  this  area  has  been  described  in  several  previous 
publications  and  symposia  presentations  (References  1-16). 

Impurity  Complexes 

Most  recently,  we  have  been  investigating  the  behavior 
of  five  transitions  noticeably  different  from  the  simpler 
transitions  discussed  earlier.  The  five  lines  are  marked  A 
through  E  in  Figure  2,  and  the  measured  energy  of  each 
transition  is  also  noted.  As  can  be  seen  in  that  figure,  the 
transitions  appear  to  have  energy  widths  similar  to  bound 
excitons  and  considerably  narrower  than  the  transitions 
involving  band  states.  It  can  also  be  observed  that  the 
energies  of  the  transitions  seem  to  indicate  rather  loosely 
bound  levels. 

Greater  insight  into  the  nature  of  the  states  involved 
in  these  transitions  can  be  gained  by  observing  their  behavior 
in  the  presence  of  external  perturbing  fields.  Our  facility 
permits  the  application  of  well  controlled  magnetic  fields  up 
to  45,000  gauss,  and  recently,  we  have  begun  exploiting 
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stress  fields  of  up  to  400  g/nini“2  which  can  be  created 
during  sample  growth  or  mounting.  Interactions  between  the 
mobile  charge  carriers  (electrons  and  holes)  and  the  magnetic 
field  causes  measurable  changes  in  the  energies  of  initial 
and  final  states  involved  in  the  transitions.  Similarly 
lattice  perturbations  due  to  applied  stresses  similarly 
create  fields  which  alter  the  energies  of  initial  and  final 
transition  states.  The  result  of  perturbations  in  the 
energies  of  the  carrier  states  is  reflected  by  changes, 
generally  encompassing  energy  shifts  and  line  splitting,  in 
the  luminescence  spectra.  The  behavior  of  the  five  specified 
lines  in  the  presence  of  magnetic  fields  is  indicated  in 
Figure  3.  The  shift  in  the  "center  of  mass"  of  transitions 
has  been  removed  in  order  to  emphasize  the  nature  of  split¬ 
ting  as  a  function  of  applied  magnetic  fields.  As  can  be 
observed  the  splitting  over  this  range  of  fields  is  linear 
(Zeeman  effect);  while  the  diamagnetic  shift  not  shown  here 
is  generally  proportional  to  higher  order  terms  (principally 
field  intensity  squared)  and  is  usually  observed  as  an 
increase  in  the  mean  transition  energy  with  applied  field. 

Our  analysis  of  the  multiplicity  of  the  lines  in  magnetic 
fields,  of  the  symmetry  of  the  splitting  patterns,  of  the 
relative  amplitudes  of  the  single  lines  in  each  pattern,  and 
of  the  magnitudes  of  the  energy  splittings  has  identified  the 
source  of  each  of  these  observations  as  a  transition  between 
a  state  involving  an  unpaired  el ectroit  wi th  one  involving  an 
unpaired  hole.  Available  data  on  complexes  A,  C,  and  E  is 
relatively  complete.  Complexes  B  and  D,  however,  have  been 
observed  less  frequently  and  neither  is  as  intense  as  the 
others,  hence  conclusions  about  these  lines  are  more  tentative. 
Calculation  of  the  effective  Lande  splitting  factors  for 
these  complexes  has  determined  that  in  all  cases  the  electrons 
involved  behave  comparably  to  those  observed  with  the  normal 
bound  exciton  transitions  which  have  been  well  studied  (gg 
=  -0.5).  For  complexes  A,  E,  and  probably  B,  effective  Lande 
factors  determined  for  the  associated  holes  are  also  similar 
to  those  measured  for  other  excitonic  systems,  k  =  +.5.  For 
complexes  C  and  probably  D,  however,  a  considerably  larger 
hole  splitting  factor  is  required  to  explain  the  observed 
splitting;  the  value  for  complex  C  being  k  =  +2.0. 

The  behavior  of  these  complexes  in  the  presence  of  an 
applied  stress  is  shown  in  Figure  4.  Also  included  in 
Figure  4,  for  reference,  is  the  observed  stress  field  behavior 
of  the  shallow  impurity  bound  excitons.  Applied  stress 
generally  results  in  a  strain  splitting  of  the  valence  band 
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(hole)  states.  In  the  exciton  systems  for  both  acceptors  and 
donors,  the  magnitude  of  this  splitting  is  seen  to  be  rela¬ 
tively  large.  For  the  five  complexes  being  investigated, 
such  splitting  is  non-existent  in  C,  D  and  E  (behavior  of 
complex  D  is  less  certain)  and  only  about  half  as  large  in  A 
and  B.  This  observation  requires  explanations  which  do  not 
involve  hole  states  such  as  would  be  present  in  bound  exciton 
or  shallow  acceptor  impurities. 

In  our  analysis  of  the  body  of  data  available  concerning 
the  five  transitions  of  interest,  we  have  applied  a  process 
of  elimination  to  arrive  at  the  model  we  feel  most  likely 
accounts  for  our  observations.  Simple  impuri ty-to-band  tran¬ 
sitions  have  been  eliminated  due  to  the  narrow  energy  widths 
of  the  observed  lines.  Impurity  bound  excitons  have  been 
eliminated  due  to  wel 1 -documented  identification  of  such 
transitions  for  all  available  impurity  species  (References  16 
and  17,  for  example)  and  by  the  unusual  behavior  of  the  lines 
in  strain  field  environments.  The  next  simplest  complexes 
which  can  be  imagined  in  a  semiconductor  are  indicated  schema¬ 
tically  in  Figure  5.  Of  these  three  possibilities  the  neutral 
donor-neutral  acceptor  complex  (Figure  5c)  is  considered 
least  likely  since  more  complicated  splitting  behavior  would 
be  expected  in  a  strain  field  and  since  splitting  in  zero 
field  environments  due  to  exchange  interactions  between 
initial  states  might  reasonably  be  expected.  The  double 
acceptor-neutral  donor  complex  (Figure  5a)  is  also  less 
likely  since  the  observed  magnetic  field  behavior  is  much 
less  complex  than  what  would  reasonably  be  expected  due  to 
the  complicated  initial  state  of  the  transition. 

The  double  donor-neutral  acceptor  complex  (Figure  5b)  can 
be  employed  to  explain  all  aspects  of  the  data  we  have  col¬ 
lected,  with  the  single  stipulation  that  the  acceptor  in¬ 
volved  is  a  relatively  tightly  bound  impurity.  Our  observa¬ 
tion  is  further  supported  by  recent  theoretical  work  reported 
by  Dr.  Van  Vechten  at  IBM  (Reference  18)  in  which  he  has 
predicted  that  complexes  involving  anti-site  defects  should 
be  relatively  common  due  to  the  growth  conditions  for  GaAs. 
Anti-site  defects  are  crystalline  imperfections  in  which  one 
specie  of  the  compound  material  occupies  a  lattice  position 
in  which  the  other  atomic  specie  should  normally  be  found. 

In  gallium  arsenide,  arsenic  atoms  sitting  on  gallium  sites 
constitute  anti-site  defects  which  would  behave  like  double 
donor  impurities.  Our  observation  of  photoluminescence 
transitions  which  support  this  prediction  are  the  first 
experimental  verification  of  their  existence  in  GaAs. 
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Implications  for  Device  Development 

The  development  of  non-destructive  characterization 
techniques  and  their  application  to  the  understanding  of 
materials  for  advanced  device  applications  are  important 
contributions  to  the  technology  base  from  which  future  sys¬ 
tems  will  evolve.  Identification  and  understanding  of  defect 
states  in  gallium  arsenide  provide  important  information, 
affecting  decisions  concerning  the  directions  of  device 
development  efforts.  As  control  of  unwanted  impurities  in 
gallium  arsenide  becomes  increasingly  well-established, 
ultimate  performance  characteristics  for  GaAs  devices  will 
become  more  dependent  upon  inherent  defects  in  the  material. 
The  understanding  of  such  defects  which  evolves  from  the 
present  work  will  provide  the  basis  for  extension  of  future 
technology  to  meet  the  ever  increasing  demands  confronting 
Air  Force  electronic  systems. 
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COMPLEXES 

COMPOSITE  PHOTOLUMINESCENCE  SPECTRUM 


FIGURE  2 


GoAs  COMPLEXES.  LINEAR  MAGNETIC  FIELD  SPLITTING  (Diamagnetic  shift  removed) 
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FIGURE  4 


Double  Acceptor  Double  Donor  Donor 
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Microcircuic  Analysis  Techniques  Using  Field _ Effecc  Liquid  Crystal 

Abscract 

Recenc  advances  of  microeleccron ics  technology  have  resulted 
in  the  availability  of  extremely  sophisticated  and  complex 
devices.  Large  scale  integration  (LSI)  circuits  achieve  a 
higher  functional  density  with  better  performance  and 
potentially  much  lower  costs.  The  reliability  of  these 
devices  is  a  key  factor  in  their  widespread  implementation  in 
military  systems.  However,  the  complexity  of  these  parts 
makes  reliability  analysis  of  these  devices  very  difficult. 
New  techniques  are  required  for  product  evaluation, 
electrical  characterization,  failure  analysis  and  reliability 
assurance. 

The  new  analysis  techniques  described  in  this  paper 
promise  to  greatly  reduce  the  time  and  effort  required  to 
pinpoint  the  cause  of  failure  in  LSI  devices  which  have 
failed  in  application  or  in  laboratory  testing.  These 
techniques  use  liquid  crystal  display  technology  to  allow  the 
analyst  to  observe  with  a  standard  optical  microscope  the 
flov7  of  data  in  a  complex  circuit  as  it  operates. 

One  method  produces  an  optical  display  of  selected  parts 
of  a  complex  circuit  by  controlling  the  electric  fields 
produced  in  a  layer  of  liquid  crystal  on  the  circuit.  This 
is  accomplished  by  special  modification  of  the  external 
control  inputs  to  the  circuit  as  the  circuit  is  operated. 
The  method  is  useful  in  design  analysis  and  chip  mapping 
where  functional  blocks  and  individual  cells  must  be 
physically  located  in  a  large  circuit.  Another  methqd  uses 
this  modified  timing  and  an  externally  applied,  synchronous 
electric  field  to  produce  a  static  display  of  dynamic  logic 
levels  present  at  internal  circuit  nodes  as  the  device  is 
operated  at  a  high  clock  rate.  This  method  is  useful  for 
nondestructive  testing  of  internal  circuit  paths  and  for 
analysis  of  low  voltage  technologies. 

Many  microcircuit  manufacturers  and  major  users  who 
maintain  microcircuit  reliability  analysis  laboratories  are 
now  implementing  and  evaluating  these  analysis  techniques. 
They  are  being  used  at  RADC  in  the  reliability  assessment  of 
new  and  advanced  microcircuit  tehnologies. 
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Introduction 


The  availability  of  .high  reliability  nicrocircui ts  to 
the  military  electronics  systems  designer  is  assured  through 
a  system  of  .procurement  .specification  docum.ents  and  by 
rigorous  reliability  evaluations  of  new  technologies  ns  ’;they 
emerge.  .Innbyations  in  .integrated  circuit  manufacturing 
methods  and  improvements  in  process  controls  "have  /resulted  dn 
microcircuits  iwhich  contain  literally  thousands  of  active 
elements  on  silicon  chips  measuring  less  than  a  centimeter  on 
a  side.  A  key  task  in  the  reliability  assessment  of  these 
complex  large  scale  .integration  (LSI)  devices  is  the 
development  of  techniques  which  aTldw  in-depth  /design 
analysis  and  accurate  failure  analysis  of  parts  vv;hich  lhave 
failed  during  laboratory  testing  or  in  ifield  use.  The 
identification  of  tontr ibut  ing  failure  mechanisms  .is 
necessary  in  order  to  recommend  device  changes,  -process 
changes  or  special  screen  tests  which  will  eliminate  the 
specific  /problem  in  future  devices. 

/Highly  cbmplex  devices  .pose  a  major  /problem  to  the 
analyst  .in  failure  verification  and  design  analysis. 
Established  techniques  rely  on  external  parametric 
measurements  and  internal  electrical  measurements  made  in 
conjunction  with  an  analysis  of  the  schematic  diagram  and 
physical  layout  of  the  ■circuit.  Circuit  diagrams  have 
generally  been  available  for  small  scale  inegration  chips, 
and  can  be  easily  verified  by  optical  inspection  of  -the 
silicon  chip.  The  number  of  active  elements  for  these  .types 
of  circuits  is  typically  less  than  a  couple  .hundred. 
Electrical  die  probing  with  micromanipulators  and  -scanning 
electron  microscope  or  -optical  inspection  techniques  are  used 
routinely  to  isolate  the  cause  of  failure  in  such  circuits. 
However.,  these  techniques  alone  are  not  adequate  tor  practical 
for  circuits  ranging  lup  to  '5,0  00  active  /eiements.  VCircuit 
diagrams  for  these  complex  'Tcircuits  are  not  generally 
available  (or  up  to  date)  and  even  the  correlation  of  circuit 
function  to  /physical  location  /on  these  devices  becomes  a 
formidable  task.  Such  information  as  v/ell  as  /other  details 
of  circuit  implementation  :may  not  be  (available  to  /the 
analyst.  Under  these  conditions,  the  only  practical  approach 
to  a  meaningful  analysis  is  with  a  /combination  of  software 
fault  isolation  techniques  and  circuit  component  level 
analysis  carried  out  in  a  •relatively  /small  .portion  of  rt’he 
circuit.  First  the  malfunction  must  be  isoiated  to  a 
particular  part  of  the  -circuit  and  then  that  part  /must  be 
analyzed  in  detail  using  techniques  similar  to  those  used  now 
as  well  as  new,  techniques  which  are  especially  suited  for 
nondestructive  testing  of  complex,  densely  packed  LSI 
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c ir cults 


The  analysis  techniques  described  in  this  paper 
represent  new  improvements  to  a  novel  analysis  technique 
originally  described  by  D.  Channin^  for  observing  the 
operation  of  integrated  circuits  using  field  effect  liquid 
crystals  and  an  optical  microscope.  Chahnin's  technique  is 
capable  of  detecting  localized  heating,  logical  "stuck  at" 
faults,  and  for  observing  data  propagating  through  operating 
circuits  at  low  frequencies  (approximately  5  Hz). 

The  present  work  was  undertaken  to  apply  the  basic 
technique  to  the  problem  of  determining  ,  the  logic  levels 
present  at  the  internal  circuit  nodes  of  a  complex  LSI  device 
during  selected  digital  states  as  the  device  under  test  (DUT) 
is  operated  at  a  high  clock  rate  in  some  predetermined  test 
sequence.  The  resulting  techniques  take  advantage  of  the  RMS 
threshold  voltage  dependence  of  the  el ectro“o pt ic  effect  in 
the  liquid  crystal  to  produce  a  display  which  is  similar  to 
that  obtained  using  scanning  electron  microscope  (SEM) 
techniques  with  beam  strobing  or  synchronous  DUT  b ias . 2 , 3 , 4 

These  techniques  are  capable  of  gate  and  transistor 
level  analysis  and  can  aid  in  fault  isolation  and  failure 
analysis.  They  are  relatively  simple  and  inexpensive  to 
implement  and  have  been  used  as  an  alternative  to  destructive 
mechanical  die  probing  and  as  a  complementary  technique  to 
SEM  voltage  contrast  methods  for  observing  device  operation 
and  determining  the  physical  location  of  various  circuit 
functional  blocks.  They  do  not  require  complex  raster 
scanning  or  sample  and  hold  elecrtonics  as  do  some  SEM  and 
laser  photoscan  analysis  techniques. 

A  description  of  chip  display  cell  preparation  is  given 
and  methods  for  generating  the  required  DUT  control  are 
discussed.  A  versatile  microcomputer  based  system  for  DUT 
control  has  been  implemented  and  used  for  these  LCD 
investigations  as  well  as  synchronous  bias  SEM  analyses. 5 
Typical  results  are  shown  which  illustrate  applications 
including  visual  determination  of  metallization 
discontinuities>  observation  of  data  associated  v/ith  selected 
digital  states  and  memory  chip  mapping. 

Display  Cell _ Prepara  tion 


These  LCD  techniques  employ  an  el 
effect^  which  depends  on  localized  changes 


ec  tro~o 
in  t  he 


ptic  display 
re  fract ive 


index  of  a  thin  layer  of  liquid  crystals.  The  integrated 
circuit  display  cell  described  by  Channin^  uses  the  DUT 
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surface  with  scratch  protection  layer  intact  and  a  thin 
glass  coverplate  to  confine  a  thin  layer  of  nematic  liquid 
crystal  material  such  as  MBBA^  as  shown  in  the 
cross-section  in  Figure  1.  The  required  vertical  molecular 
orientation  can  be  promoted  in  the  liquid  crystal  by  first 
treating  the  surfaces  with  a  surfactant  solution. 
Excellent  results  have  also  been  obtained  using  a  surface 
treatment  described  by  Janning  for  the  coverplate  which 
uses  a  thin,  transparent  Au  layer  which  is  sputtered  onto 
the  coverplate  surface.  The  coverplate  must  be  cut  with  a 
v;afer  scriber  to  fit  within  the  perimeter  defined  by  the 
bonding  wires  so  that  it  does  not  interfere  with  the 
electrical  connections  to  the  device.  A  drop  of  liquid 
crystal  is  then  placed  on  the  device  surface.  The 
coverplate  is  then  floated  on  the  drop  of  liquid  crystal. 
A  typical  prepared  cell  is  shown  in  Figure  2. 


The  device  is  then  fixtured  under  a  standard 
microscope  and  viewed  with  vertical  illumination.  A 
polarizer  is  placed  in  the  light  path  and  another  is  placed 
in  the  viewing  path  and  oriented  at  90°  to  the  first.  The 
display  effect  cannot  be  observed  without  vertical 
illumination  and  properly  oriented  polarizers.  After  an 
analysis  is  completed,  the  liquid  crystal  can  be  removed  by 
cleaning  the  device  with  acetone  and  the  lecithin  coating 
can  be  removed  with  Uresolve  Plus^. 


Display  Cell  Characteristics 

The  display  cell  formed  as  above  is  a  field  effect 
display,  in  contrast  to  the  dynamic  scattering  display 
described  by  Salvo^  for  detecting  defects  in  insulating 
films.  When  the  scratch  protection  glass  is  left  intact 
the  cell  shows  no  display  effect  for  dc  voltages  as  a 
result  of  surface  charging  (with  a  time  constant  of  about 
0.1  sec)  within  the  display  cell  which  reduces  the  vertical 
component  of  electric  field  in  the  liquid  crystal  below  the 
threshhold  level.  (V7hen  the  glass  is  removed,  dc  voltages 
can  be  observed).  VJhen  viewed  through  crossed  polarizers, 
the  device  appears  dark  if  no  electrical  bias  is  applied, 
as  shown  in  Figure  3.  Note  that  dust  particles  in  the  cell 
cause  bright  spots.  These  spots  can  be  distinguished  from 
active  circuit  elements  since  they  do  not  change  with  the 
removal  and  application  of  bias.  When  bias  is  applied 
and  the  device  is  operated  there  is  a  bright  display  over 
those  conductors  carrying  sufficiently  high  RMS  voltages 
{ 4 . 5v  RMS).  The  display  is  continuous  for  frequencies 
above  500  llz  to  greater  than  10  MHz.  Figure  4  shows  the 
display  for  a  14  stage  binary  counter  (CMOS,  Vj^j^  =  +io  V) 
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clocked  at  a  10  MHz  race.  The  14  D-cype  flip-flop  stages 
are  all  toggling  with  a  50%  duty  cycle  at  different 
frequencies.  The  voltage  waveforms  carried  by  the  signal 
conductors  in  each  of  the  stages  have  the  same  RMS  value 
(7.07  V)  and  are  all  displayed. 

Device  Control  For  A  Selective  Display 

The  display  produced  by  the  counter  circuit,  shown  in 
Figure  4  does  not  give  information  about  the  logic  levels 
present  on  various  conductor  lines  during  a  g iven  • d ig ica 1 
state.  The  RMS.  voltage  values  are  the  result  of  an 
effective  integration  of  the  logic  levels  present  during 
the  period  of  a  complete  cycle  through  all  of  the  16384 
individual  digital  states.  The  display  is  continuous 
because  the  complete  cycle  is  repeated  every  1.6+  ms  for 
the  10  MHz  clock  rate  used  in  this  example.  It  is 
desirable/  however/  to  produce  a  display  which  does  give 
information  about  the  logic  levels  present  during 
individual  digital  states  in  the  test  cycle  and  also  about 
the  physical  location  of  various  circuit  functional  blocks/ 
such  as  single  flip-flops  or  clock  bus  lines. 

A  display  y^hich  does  provide  the  desired  information 
is  shown  in  Figure  5.  The  circuit  is  operated  in  quite  a 
different  manner  to  produce  the  selective  display.  The 
timing  of  state  transitions  in  the  test  cycle  has  been 
modified  such  that  the  circuit  "pauses"  in  two  preselected 
states  for  a  time  much  greater  than  the  duration  of  any  of 
the  other  states  in  the  cycle.  This  modified  timing 
effectively  increases  the  duty  cycles  of  two  states  in  the 
test  cycle  and  reduces  the  duty  cycles  of  all  the  remaining 
states  in  the  test  cycle.  As  a  result/  the  RMS  voltage 
values  on  the  conductors  on  the  chip  become  largely  a 
function  of  the  logic  levels  present  only  during  the  tv/o 
preselected  states.  The  display  contains  a  bright  response 
only  over  those  conductors  which  assume  different  logic 
levels  in  the  two  chosen  states.  This  special  test  cycle 
timing  is  illustrated  in  Figure  6.  In  this  manner  a 
"reference"  state  (e.g.  binary  0000)  and  a  "selected"  state 
(e.g.  binary  1000)  can  be  designated/  and  the  logic  levels 
present  during  any  single  selected  state  in  the  test  cycle 
can  be  directly  compared  to  those  present  during  the 
reference  state.  The  display  in  Figure  5  allows  the 
analyst  to  easily  determine  the  location  of  the  fourth 
flip-flop  cell  and  also  to  verify  the  correct  logical 
content  of  the  flip-flop  during  the  selected  state. 
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It  IS  interesting  to  note  that  information  very  similar  to 
SEM  methods  is  obtained  without  a  busy  background  of  dc 
level  information  unrelated  to  the  feaiiur6s  of  interest. 
This  can  simplify  interpretation  of  the  display  for  lov/ 
magnification  searching  such  as  that  required  for  chip 
mapping  of  a  memory  array  or  for  circuit  conductor  tracing 
when  the  conductor  path  is  distributed  tihrbughout  large  and 
widely  separated  areas  on  the  chip. 

The  technique  described  above  <lan  be  used  to 
investigate  a  variety  of  circuit  functions;  shift 
registers/  control  logic  and  any  circuitry  which  can  be 
controlled  in  the  required  manner.  The  basic  idea  is  to 
accentuate  certain  information  by  causing  selected 
conductors  to  change  logic  levels  at  a  near  50%  duty  cycle 
and  to  suppress  other  information  by  decreasing  the  duty 
cycle  (and  RMS  value)  of  the  logic  levels  associated  with 
those  states  not  to  be  displayed. 

The  control  signals  for  producing  the  pauses  in  the 
periodic  test  cycle  can  be  generated  in  a  variety  of  v\'ays. 
Analog  delay  circuits  (e.g.  one  shots)  or  digital  delays 
(e.g.  counters)  can  be  used  to  gate  a  clock  input  for  many 
circuit  functions.  These  approaches  have  the  disadvantage 
that  they  require  hardwired  circuits  which  are  unique  for 
each  circuit  function  and  number  of  states  in  the  test 
cycle  and  are  not  easily  implemented  when  more  than  one 
control  line  is  required.  A  more  versatile  approach  which 
incorporates  a  small/  one  board  microprocessor  system  to 
generate  DUT  control  signals.  The  OUT  timing  diagram 
signals  are  converted  to  a  sequence  of  bytes  which  are 
stored  in  the  microcomputer's  memory  and  outputted  under 
software  control.  The  appropriate  pauses  are  thus  easily 
programmed  into  the  cycle  and  can  be  changed  by  the 
programmer  or  automatically  by  the  program.  Applications 
which  use  this  additional  capability  will  now  be  discussed. 
(The  microcomputer  used  for  this  work  was  based  on  the  1{j02 
CMOS  device/  with  only  512  bytes  of  program  storage.  Total 
cost  of  the  system  was  less  than  $400.00). 

Applications  and _ Results 

A  simple  but  very  usef\il  application  of  the  LCD 
technigue  is  to  solve  the  problem,  of  verifying  continuity 
of  conductor  lines  at  internal  nodes  of  small  geometry  LSI 
devices  where  mechanical  die  probing  Would  be  difficult  at 
best.  Figure  7  shows  a  display  which  identifies  the  exact 
location  of  electrical  opens  which  were  the  result  of 
corrosion.  This  device  was  a  field  failure.  The  circuit  is 
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a  1024  bic  PMOS  static, shift  register.  The  opens  occurred 
on  clock  bus  lines  and  resulted  in  catastrophic  failure. 
The  portions  of  the  conductor  lines  which  are  dark  in 
Figure  7  are  electrically  not  connected  to  the  IKHz  input 
clock  or  the  other  two  clock  phases  which  are  internally 
generated.  Destructive  mechanical  die  probing  was  not 
necessary  in  this  case  to  identify  the  fail  sites. 

A  task  v/hich  occurs  frequently  in  the  analysis  of  LSI 
devices  is  chip  mapping.  For  the  Case  of  a  memory  circuit, 
the  binary  address  corresponding  to  each  bit  stored  must  be 
correlated  to  a  physical  memory  cell  location  on  the  chip. 
This  usually  requires  that  the  address  decode  circuitry  be 
completely  traced  to  identify  ceil  locations.  Lov/ 
magnification  searching  with  optical  or  SEM  techniques  can 
be  tedious  and  time  consuming.  The  task  can  be  simplified 
considerably  by  using  the  LCD  technique  with  the  DUT  under 
automatic  control.  Figure  8  illustrates  the  static  display 
obtained  by  writing  an  alternate  pattern  of  1/0  into  two 
bit  locations;  the  binary  ”0"  address  as  a  reference  and 
one  other  selected  cell  which  is  then  easily  identified. 
The  memory  array  is  easily  mapped  as  the  DUT  controller  is 
sequenced  from  one  selected  address  to  the  next.  The  device 
shown  is  a  CMOS  256  bit  memory.  Similar  results  have  been 
obtained  for  an  UV  erasable  ROM  (  2708,  8K  bit)  and  other- 
memory  arrays. 

Recent  efforts  have  been  directed  at  applying  these 
LCD  techniques  to  the  analysis  of  a  microprocessor  device. 
The  device  selected  v/as  the  1802  device  v/hich  is  a  CMOS 
circuit  which  can  be  operated  on  a  single  supply  and  could 
be  easily  tested  in  the  same  microcomputer  system  which  was 
used  for  controlling  other  devices  analyzed  with  these 
techniques.  Much  information  can  be  generated  by  simply 
observing  the  LCD  display  as  the  microprocessor  executes 
test  programs  v/hich  execute  certain  functions  on  the  chip. 
However,  it  is  desirable  to  inspect  the  data  flow  in  the 
circuits  throughout  the  chip  during  the  intermediate  steps 
in  the  machine  cycles  which  time  the  execution  of  each 
instruction  in  the  test  software.  This  .required  the  use  of 
a  second  microcomputer  system  to  synthesize  the  master 
clock  and  control  signal  inputs  for  the  microcomputer 
containing  the  device  being  tested  with  a  LCD  display  cell. 
Those  control  signals  had  to  be  specially  timed  to  allow 
the  reference  state  to  selected  state  comparison  v^/hich  was 
used  previously  for  smaller  circuits. 

Videotape  recording  is  used  to  capture  the  displays 
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resulting  from  the  execution  of  a  test  program  which 
exercises  certain  portions  of  the  device  with  different 
data  and  instructions.  The  controlling  microcomputer  can 
be  programmed  to  dwell  at  one  particular  instruction  in  the 
test  program  or  to  sequence  from  one  to  the  next/  so  that 
one  state  can  be  inspected  for  a  long  time  or  a  long 
series  of  states  can  be  recorded. 

The  following  motion  picture  shows  the  display 
obtained  for  some  typical  test  programs  v;hich  exercise  the 
microprocessor's  machine  cycle  counter  and  decoder,  a 
program  counter  register  (RO),  and  finally  the  DATA  (or 
accumulator)  register. 

(  Kotion  picture  running  time  approximately  5  minutes). 

Detailed  comparisons  can  be  made  of  tv/o  states  of  the 
device  and  correlated  to  the  execution  of  various 
instructions  and  physical  locations  on  the  chip.  For 
example,  a  photomicrograph  may  be  produced  by  making  an 
exposure  filtered  in  green  and  an  exposure  filtered  in  red, 
each  with  a  different  selected  state.  Those  conductors 
which  assume  different  logic  levels  during  both  of  these 
selected  states  than  the  logic  levels  assumed  during  the 
reference  state  show  in  a  third  color.  Thus  data  which 
changes  from  state  to  state  as  instruction  execution 
proceeds  can  be  highlighted.  Those  conductors  which  do  not 
change  from  their  initial  values  in  the  reference  state  do 
not  light  and  thus  do  not  interfere  with  the  interpretation 
of  the  active  portions  of  the  display. 

These  techniques  allow  the  analyst  to  correlate 
circuit  function  to  physical  location  on  the  chip  without 
the  aid  of  a  complete  chip  map  and  circuit  schematic.  The 
contents  of  internal  registers  and  bus  lines  can  be 
inspected  for  different  states  of  the  circuit.  Improper 
operation  can  be  detected  by  comparison  between  the 
displays  for  different  states  or  by  compairision  to  the 
display  of  a  similar  good  device.  . 

Biased  Coverplate  Method 

Preliminary  results  have  been  obtained  which  have 
proven  the  feasibility  of  depositing  a  conductive  layer  of 
Au  on  the  lower  surface  of  •,  the  coverplate  and  making 
electrical  contact  to  it  to  form  a  field  plate.  The  field 
plate  voltage  can  then  be  used  to  directly  control  the 
digital  state  selected  for  display.  The  advantages  of  this 
approach  are  that  the  weaker  electric  fields  available  on 
devices  with  small  voltage  swings  can  be  aided  by  pulsing 


the  field  place.  For  example/  a  simple  MAND  gate  device 
may  be  biased  with  a  4v  square  wave  on  one  input.  The 
conductive  coverplate  is  biased  with  a  4v  (Square  wave  IBO 
degrees  out  of  phase  with  the  input  waveform.  In  the 
resulting  display  those  conductors  which  carry  logic  levels 
out  of  phase  with  the  fieldplace  are  imaged.  This  method 
also  allows  display  of  only  the  logic  0  conductors  or  only 
the  logic  1  conductors  for  the  selected, state .  The  phase  of 
the  field  plate  waveform  determines  which  logic  level  is 
displayed.  Experiments  to  date  have  used  a  conductive  epoxy 
to  establish  contact  to  the  field  place  on  the  lov/er 
surface  of  the  coverplate.  The  epoxy  is  attacked  by  the 
liquid  crystal  and  contact  is  eventually  lost.  Presently, 
experiments  are  being  made  with  a  fired  on  Au  layer  which 
provides  a  metal  path  around  the  edge  of  the  coverplate  to 
the  top  surface  which  can  then  be  contacted  from  above  with 
a  micromanipulacor  probe.  This  biased  coverplate  method 
extends  the  technique  to  the  lower  voltage  technologies 
such  as  NMOS  and  TTL. 

Conclusions 

The  LCD  analysis  techniques  described  give  the  analyst 
another  tool  which  can  be  used  routinely  in  addition  to  the 
many  available  techniques  already  being  applied  to  the 
problems  encountered  in  the  analysis  of  complex  digital 
devices.  The  materials  required  are  commercially  available 
and  the  equipment  needed  is  likely  to  be  present  in  the 
analysis  laboratory  or  available  at  relatively  low  cost. 
The  techniques  are  nondestructive  and  the  analyst  can 
easily  remove  the  LCD  cell  after  use  and  proceed  to  other 
steps  such  as  die  probing  or  SEM  analysis. 

The  results  which  are  obtained  with  these  techniques 
are  similar  to  those  for  SEM  and  laser  photoexcitation 
techniques^®  which  require  beam  modulation,  T.V.  raster 
scanning,  and  complex  display  electronics.  Those  latter 
techniques  may  represent  a  prohibitively  high  initial 
investment  cost  to  the  smaller  laboratory  and  also  require 
specially  trained  personnel  to  operate  and  maintain  the 
required  equipment. 

The  LCD  display  can  be  made  to  respond  only  to  ac 
voltages  when  the  chip  overcoat  glassivatioh  is  left 
intact,  ih  contrast  to  SEM  techniques  which  also  image  dc 
voltages.  This  feature  is  used  to  simplify  the 
interpretation  of  the  LCD  display;  especially  when 
searching  for  small  features  at  low  magnification.  The 
actual  active  display  areas  are  hiqh  contrast,  light  on 
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dark  and  are  confined  to  the  immediate  area  above  conductor 
lines  on  the  surface  of  the  circuit.  This  type  of  display 
is  quite  different  from  that  obtained  v/ith  photoexcitation 
and  electron  beam  induced  current  (EBIC)  techniques!!  which 
are  used  primarily  to  image  p-n  junctions.  The  LCD 
technique  does  not  give  structural  information  and  cannot 
be  used  to  identify  the  location  of  an  inversion  region, 
for  example,  as  can  some  of  the  previously  mentioned 
techniques.  However>  the  display  is  easily  correlated  to 
the  optical  image  of  the  circuit  by  simply  rotating 
(uncrossing)  one  of  the  polarizers  in  the  optical  viewing 
system. 

Finally,  the  implementation  and  interpret 
these  LCD  techniques  are  easily  within  the  capabil 
constraints  of  the  modestly  equipped  microcircuit 
facility. 
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Figure  1 

Cros s-s ecc ion  of  a  display  cell. 
The  horizoncai  allighmenc  in  the 
liquid  crystal  above  the  third 
conductor  on  the  right  is  in 
response  to  fringing  fields  from 
the  underlying  conductor. 
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Display  . cell  prepared  on  a  CMOS 
binary  counter.  (4020). 


Selective  display  of  seventh  stage  b 
comparing  states  "0"  and  "64". 
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Figure  6 

Control  cydle  timing  for  a  selective 
display  of  fourth  stage  by  comparing 
states  "0"  and  "8". 


Figure  B 

Reference  memory  cell  location 
is  at  lower  right.  Selected  c 
"81"  is  at  upper  left. 


Electrical  opens  on  clock  lines 
this  1024  bit  PMOS  static  shift 
register  are  identified  by  the 
'discontinuous  displays. 


I 


of 


Biocrraphical  Sketch 


Daniel  J.  Burns  was  born  in  Grand  Rapids /  Michigan  on 
April  7/  1950.  He  received  a  BEE  degree  from  the  University 
of  Detroit»  Detroit/  MI  in  1973  and  the  MSEE  degree  from 
Syracuse  University,  Syracuse,  NY  in  1978. 

Since  1973  he  has  been  with  the  Reliability  Branch  at 
Rome  Air  Development  Center,  Griffiss  AFB,  NY  where  his 
primary  work  involves  microcircuit  analysis  in  support  of  DOD 
systems.  He  has  also  worked  on  reliability  evaluations  of 
advanced  microcircuit  technologies  and  on  specifications  for 
microcircuit  electrical  overstress  tolerance.  He  is 
currently  engaged  in  development  of  techniques  for  analysis 
of  complex  devices  using  liquid  crystals. 


84 


SURFACE  ACOUSTIC  WAVE  FREQUENCY  SYNTHESIZER  FOR  JTIDS 


BY 


ALAN  J.  BUDREAU,  ANDREW  J.  SLOBODNIK,  JR.,  and  PAUL  H.  CARR 


Electromagnetic  Sciences  Division 


Rome  Air  Development  Center 
Depyty  for  Electronic  Technology 
Hanscom  AFB,  Massachusetts  01731 


85 


Surface  Acoustic  Wave  Frequency  Synthesizer  For  JTIDS 


Abstract 


Fast-hopping  frequency  synthesizers  are  essential  to 
a  variety  of  communication  and  data-link  applications.  A 
frequency  synthesizer  built  as  a  feasibility  demonstration 
for  the  Joint  Tactical  Information  Distribution  System 
(JTIDS)  is  described.  Surface  Acoustic  Wave  (SAW)  tech¬ 
nology  was  utilized  to  yield  a  design  amenable  to  a  com¬ 
pact,  reliable,  low  cost  final  product.  Although  any  of 
the  standard  approaches  to  frequency  synthesis,  i.e., 
phase-locked-loop  oscillators,  digital  synthesis,  and 
direct  synthesis,  may  benefit  from  SAW  components  to 
reduce  the  size  and  cost,  the  direct  synthesizer  is  the 
most  appropriate  for  achieving  fast  frequency-hopping. 

The  JTIDS  terminals  require  sub-microsecond  hopping. 

A  mixing  approach,  in  which  the  outputs  of  three  direct 
synthesizers  are  mixed,  was  chosen  to  reduce  the  number 
of  filters  required.  Advanced  state-of-the-art  UHF  SAW 
filters  were  developed  and  combined  with  low  loss  multi¬ 
plexing  techniques  and  switches  into  high  performance 
miniature  SAW  filterbanks.  Design,  fabrication  and 
testing  are  described.  Data  from  the  prototype  frequency 
synthesizer  is  presented. 
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Introduction 

The  purpose  of  this  work  was  to  develop  new  techno¬ 
logy  for  a  low  cost,  compact  frequency  synthesizer  for 
integrated  terminals  for  the  Joint  Tactical  Information 
Distribution  System  (JTIDS)  and  Tactical  Air  Navigation 
System  (TACAN) .  Such  termirial^  require  at  least  2  simul¬ 
taneous  independent  tones  chosen  from  among  252  frequen¬ 
cies.  This  is  particularly  applicable  to  fighter  air¬ 
craft,  where  size,  weight,  and  cost  savings  are  highly 
significant.  Based  on  a  total  of  5,000  units,  savings  on 
the  order  of  $20  million  (unit  cost  reduced  from  $6000  to 
$2000)  could  be  realized.  This  cost  reduction,  along  with 
achieving  a  size  of  60  cubic  inches  or  less,  are  the 
expected  results  of  a  recently  awarded  Manufacturing  Tech¬ 
nology  (ManTech)  contract  based  on  the  present  work. 

Design  Approach 

Fractional  microsecond  frequency-hopping  is  required 
for  JTIDS.  Frequency  synthesizers  based  on  retuning  phase- 
locked  loops  are  inherently  much  too  slow  for  this  require¬ 
ment.  The  rapid  hopping  must  be  obtained  without  sacri¬ 
ficing  quality,  i.e.,  the  output  tones  must  be  locked  to 
the  system  clock  and  of  low  noise  and  spurious  content.  At 
the  same  time,  the  intended  airborne  and  manpack  uses  dic¬ 
tate  that  size,  weight,  and  power  consumption  be  minimized. 
The  unit  to  be  described  here  is  a  prototype  breadboard 
designed  to  demonstrate  feasibility  by  synthesizing  219  of 
the  tones.  The  ManTech  contract  will  improve  the  synthe¬ 
sizer  to  achieve  2  simultaneous  outputs,  all  252  tones,  and 
a  compact,  readily  manufacturable  unit. 

The  concept  of  the  surface  acoustic  wave  (SAW)  direct 
synthesizer  is  illustrated  in  Figure  1.  A  comb  generator, 
or  source  of  multiple  frequencies  as  indicated  in  the  spec¬ 
trum  analyzer  display  in  the  upper  left-hand  corner,  is 
locked  to  a  stable  clock.  The  stability  and  accuracy  of 
the  tones  is  that  of  the  clock.  The  comb  spectrum  is  simul¬ 
taneously  fed  to  the  9  channels  of  the  SAW  filterbank,  which 
sorts  them  according  to  frequency,  so  that  each  input  to  the 
9X1  switch  carries  a  single  but  different  tone.  The 
desired  tone  is  then  selected  by  the  9  X  1  switch  as  indi¬ 
cated. 

The  basic  direct  synthesizer  architecture  of  Figure  1 
could  be  realized  with  any  type  of  filter.  It  is  this 
design,  in  which  the  tones  present  at  the  output  of  each 
filter  are  continuous  waves  (CW) ,  and  therefore  capable  of 
being  switched  as  fast  as  the  switches  themselves. 
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that  yields  the  high  inherent  speed  of  frequency-hopping. 
However,  most  filter  technologies  would  lead  to  a  bulky 
filterbank.  SAW  filterbanks,  on  the  other  hand,  are  com¬ 
pact,  as  well  as  being  capable  of  manufacture  by  a  single 
photolithographic  step.  This  useful  combination  of  the 
direct  synthesizer  approach  with  SAW  technology  was  the 
subject  of  Patent  no.  3,891,929  awarded  to  Carr  and  Budreau 
in  1975. 

The  compactness  of  SAW  devices  derives  from  the  velo¬ 
city  of  the  surface  wave  being  typically  1/100,000  that 
of  the  speed  of  light.  Dimensions  of  resonant  structures 
scale  in  the  same  ratio.  Thus,  at  300  MHz,  the  wavelength 
of  an  electromagnetic  wave  is  1  meter;  the  corresponding 
SAW  has  a  wavelength  of  0.01  millimeter.  Nine  SAW  filters 
may  be  seen  as  the  dark  rectangle  on  the  left  end  of  the 
package  of  Figure  2.  Even  though  they  were  not  packaged 
for  minimal  size,  they  only  occupy  a  space  of  1.6  X  2.7  cm. 

Let  us  consider  the  operation  of  a  single  SAW  filter. 

It  is  shown  schematically  in  Figure  3.  The  electromagnetic 
input  signal  excites  an  electric  field  between  the  conduct¬ 
ing  grids  of  the  transducer  and  is  converted,  due  to  the 
piezoelectric  nature  of  the  substrate,  to  mechanical  strain 
of  the  surface.  This  strain  then  propagates  in  opposite 
directions  (i.e.,  in  both  directions  along  the  long  axis 
of  the  substrate).  The  waves  propagating  to  the  left  from 
the  input  transducer  are  absorbed  by  the  acoustic  absorber 
placed  there  for  just  this  purpose.  The  desired  waves  are 
those  shown  propagating  between  the  transducers,  which, 
when  they  reach  the  output  transducer,  generate  an  electri¬ 
cal  signal  across  the  load  "Ro”»  through  the  piezoelectric 
effect.  The  portion  of  the  waves  which  pass  beyond  the 
output  transducer  are  eliminated  by  the  absorber  placed  at 
that  end.  The  center  frequency  to  which  the  transducers 
are  resonant  is  determined  by  the  line  spacing  of  the  trans¬ 
ducers,  while  the  shape  of  the  frequency  response  curve  is 
determined  by  their  detailed  structure. 

Shaping  of  the  frequency  response,  or  "weighting",  may 
be  accomplished  in  a  number  of  ways.  Two  of  the  most  satis¬ 
factory  forms  of  weighting,  withdrawal  weighting  and  finger 
overlap  weighting  (apodizat ion) ,  are  shown  in  the  figure,, 
as  they  were  the  techniques  used  for  the  filters  in  the 
synthesizer. 

The  output  of  the  comb  generator ,. containing  all  of 
the  desired  tones,  must  be  distributed  to  all  9  SAW  filters. 
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Power  dividers  are  bulky  and  expensive.  At  the  RADC  Deputy 
for  Electronic  Technology,  we  have  developed  a  technique 
for  using  a  minimum  of  matching  elements , in  combination  with 
appropriate  choice  of  series  and  parallel  interconnections 
of  the  transducers  to  accomplish  the  power  division  in  a 
compact,  efficient,  yet  inexpensive  way.  This  approach  was 
awarded  Patent  no.  3,942,140  issued  to  Laker,  Carr,  and 
Budreau  in  1976,  as  well  as  a  Proceedings  of  the  lEEEpaper 
by  the  same  authors  in  May  1976.  For  the  current  synthesizer 
design,  we  incorporated  an  improved  version  of  this  multi¬ 
plexing,  to  be  discussed  in  the  SAW  and  switch  section  later 
on. 


The  requirement  for- the  generation  of  219  tones  dictated 
a  more  efficient  system  than  could  have  been  achieved  by 
directly  using  the  scheme  of  Figure  1  with  219  filters 
(252  tones  is  the  appropriate  number  for  an  integrated  JTIDS/ 
TACAN  terminal,  but  specification  changes  and  a  tight  sche¬ 
dule  dictated  the  demonstration  with  only  219).  Furthermore, 
such  a  "brute  force"  approach  would  have  had  other  problems, 
inclyding  the  wide  bandwidth  and  close  relative  spacing  of 
the  required  frequencies.  By  starting  with  three  direct 
synthesizers,  and  mixing  together  two  sets  of  9  tones  and 
one  set  of  3  tones,  we  have  achieved  a  capability  of 
9x9x3=  243  tones.  A  simplified  block  diagram  of  the 
direct  SAW  synthesizer  using  mixing  is  presented  as  Figure  4. 
The  mixer  synthesizer  consists  of  two  9-tohe  direct  SAW 
synthesizers  of  the  type  in  Figure  1,  plus  a  three  tone 
direct  synthesizer  using  three  phase-locked  loops.  For  con¬ 
venience  in  design,  fabrication,  and  packaging,  the  synthe¬ 
sizer  is  divided  into  three  sections.  The  source  section 
contains  comb  generators  covering  321-345  MHz  in  3  MHz 
increments  and  360-468  MHz  in  12  MHz  increments,  plus 
multipliers  and  phase-locked  loops  which  produce  1984,  1985, 
and  1986  MHz  tones.  In  addition,  it  has  the  3X1  switch 
which  chooses  among  the  latter  sources.  The  SAW  and  switch 
section  includes  the  two  9- tone  SAW  filterbanks,  each  of 
which  is  connected  to  a  9.X  1  switch,  just  as  in  Figure  1. 

The  output  or  mixer  section  contains  amplifiers,  filters,  two 
mixers,  and  a  frequency  doubler. 

Let  us  examine  the  hardware  required  in  more  detail,  as 
shown  in  Figure  5.  Note  that  219  tones  are  obtained  by  the 
use  of  only  18  SAW  filters  plus  the  5  frequency  sources.  '  Use 
of  the  mixing  sc'reme,  plus  frequency  doubling  to  achieve  the 
high  frequency  set  of  tones,  allows  the  filters  to  operate 
in  the  321-468  MHz  range.  Since  the  line  widths  and  spacings 
of  the  SAW  transducers  are  proportional  to  the  wavelengths 
of  the  signals  to  be  passed,  keeping  the  frequencies  below 


89 


500  ^^Hz  allows  fabrication  of  the  SAW  transducers  by  an 
extension  of  conventional  optical  photolithography. 

The  overall  system  design  was  optimized  by  choosing 
the  frequencies  to  be  mixed  in  such  a  way  as  to  minimize 
the  problems  with  spurious  signals,.  This  was  done  in 
collaboration  with  our  colleagues  W.  M.  Bridge  and 
N.  Moulding  at  MITRE  Corporation,  as  we  described  at  the 
1976  Ultrasonics  Symposium.  The  problems  of  the  spurious 
signals  produced  and  the  filtering  required  to  get  rid  of 
them  were  crucial,  as  we  were  working  to  satisfy  the  restric¬ 
tive  specifications  shown  by  the  lower  curve  of  Figure  6. 

Let  us  examine  Figure  5  more  carefully.  The  synthesizer 
consists  of  the  three  blocks  shown.  The  system  clock  (1  MHz 
in  our  current  experiments)  provides  the  reference  signal  to 
all  five  phase-locked  loop  oscillators  (PLLOs)  of  the  source 
package.  Each  of  the  sets  of  9  tones  is  used  to  generate  a 
single  signal,  as  described  in  connection  with  Figure  1. 

These  signals  are  amplified,  then  the  360-468  MHz  tone  is 
doubled  to  yield  a  720-936  MHz  tone.  The  321-345  MHz  tone 
and  the  720-936  MHz  tone  are  mixed  to  yield  a  difference  in 
the  378-615  MHz  range.  This  difference  frequency  is  again 
mixed  with  the  1984-1986  MHz  tone  to  yield  a  difference 
frequency  in  the  desired  frequency  range.  Appropriate  ampli¬ 
fication  and  limiting  are  used  to  obtain  a  leveled  output. 

The  three  blocks  will  now  be  described  in  detail. 

Source  Section 


All  the  signals  which  appear  at  the  output  of  the  fre¬ 
quency  synthesizer  arise  from  combining  tones  which  originate 
in  the  source  section;  these  are  as  stable  and  accurate  as  the 
dlock.  This  stability  and  accuracy  is  carried  through  the 
synthesizer  and  is  reflected  in  the  high  quality  of  the  out¬ 
put  signal.  On  the  other  hand,  any  noise  or  close-in  spurious 
signals  present  in  the  source  section  will  also  be  present 
in  the  output,  putting  very  tough  specifications  on  this  cir¬ 
cuitry.  A  source  package  consisting  of  the  two  comb  genera¬ 
tors  and  a  1984  MHz  source  was  obtained  under  AF  contract 
from  Zeta  Laboratories,  Santa  Clara,  CA.  The  approach  was 
to  multiply  up  from  1  MHz,  use  crystal  filters  to  select  the 
appropriate  harmonic,  and  use  this  to  drive  a  step  recovery 
diode.  The  step  recovery  diode  transforms  the  input  waveform 
into  very  narrow  .pulses  which  are  rich  in  harmonics.  RF  fil¬ 
tering  of  the  output  comb  spectrum  was  required  to  meet  the 
specifications.  The  1984  MHz  was  obtained  by  using  a  phase- 
locked  oscillator  at  99.2  MHz  locked  to  the  1  MHz  clock, 
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followed  by  a  X  20  stage.  This  was  accomplished  in  a  pack¬ 
age  of  99  cu.  in.,  which  could  be  engineered  down  to  include 
the  1985  and  1986  MHz  sources  in  the  same  package.  As  with 
other  forms  of  frequency  multiplication  the  step  recovery 
diode  approach  has  relatively  high  power  requirements. 

A  new  and  exciting  approach  for  the  source  section  is 
the  mode-locked  SAW  oscillator  (MLSO)  developed  at  United 
Technologies  Research  Center,  E.  Hartford,  CT  under  AF  con¬ 
tract.  It  was  described  at  the  1977  IEEE  Ultrasonics  Sym¬ 
posium,  and  also  in  the  final  report  under  the  contract, 

"High  Frequency  Signal  Sources  Using  Surface  Acoustic  (SAW) 
Technology,"  by  Dr.  M.  Gilden,  Report  RADC-TR-399,  1977. 

The  operation  is  shown  in  the  block  diagram.  Figure  7.  A 
recirculating  loop  including  a  SAW  delay  line,  amplifier 
(AMP  1),  and  a  non-linear  element  (AMP  2  with  low  supply 
voltage)  was  used  to  directly  produce  a  stable  comb  spec¬ 
trum.  While  a  simple  delay  loop  will  generally  "lock  on" 
to  a  single  frequency  of  oscillation,  design  of  a  delay 
path  which  allows  gain  greater  than  unity  at  the  desired 
frequencies,  plus  the  appropriate  non-linear  element  will 
yield  the  comb  spectrum.  Further,  this  multi-frequency 
oscillator  can  be  locked  to  an  external  clock.  As  this  SAW- 
based  comb  generator  has  been  shown  to  be  smaller,  lighter, 
and  to  consume  less  power  than  the  more  conventional  electro¬ 
nic  approach  described  above,  it  was  our  preferred  implemen¬ 
tation.  It  yielded  the  two  comb  spectra  plus  the  1984 
(but  not  1985  or  1986)  MHz  signals  in  a  package  of  just  under 
50  cu.  in.  This  approach  will  require  a  temperature  compen¬ 
sation  scheme  which  is  being  implemented  under  a  subsequent 
contract. 

SAW  and  Switch  Section 


In  this  section,  the  design,  fabrication  and  testing  of 
the  two  high  performance  9-channel  switchable  SAW  filterbanks 
on  temperature  compensated  quartz  will  be  described.  As 
depicted  in  Figure  5,  these  filterbanks  consist  of  9  indivi¬ 
dual  SAW  filter  channels  interconnected  to  operate  as  a 
single  input,  single  switch- selected  output,  device.  Each 
filter  consists  of  a  length  weighted  transducer  and  a  with¬ 
drawal  weighted  transducer  as  illustrated  in  Figure  3. 

The  frequency  response  of  a  single  high  frequency  filter 
shown  in  Figure  8  demonstrates  the  extremely  high  out-of-band 
rejection  of  these  SAW  filters.  In  order  to  achieve  these 
results  many  factors  must  be  optimized  and  second-order 
effects  controlled.  For  example,  diffraction  correction 
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(to  compensate  for  the  spreading  of  the  acoustic  beam)  is 
used  in  the  overlap  weighted  transducer;  velocity  correction 
(for  the  differing  acoustic  velocities  on  the  bare  substrate 
and  those  portions  covered  with  metal)  is  used  on  the  with¬ 
drawal  weighted  transducer;  tapered,  notched  substrates  are 
used  to  suppress  bulk  mode  spurious  (a  SAW  transducer  gener¬ 
ates  and  receives  volume  waves  as  well  as  SAWs) ;  and  careful 
packaging  is  used  to  minimize  electromagnetic  feedthrough. 

As  mentioned,  the  individual  transducers  must  be  inter¬ 
connected  in  order  to  achieve  the  desired  operation.  The 
goals  of  any  interconnection  or  multiplexing  technique  are 
as  follows:  1)  insertion  loss  minimization,  2)  insertion 
loss  uniformity  among  channels,  and  3)  small  size  and  low 
cost.  Results  of  our  studies  have  shown  that  the  proper 
choice  of  a  series -parallel  interconnection  scheme  can 
achieve  each  of  these  goals. 

In  particular,  input  multiplexing  is  achieved  here  by 
the  parallel  interconnection  of  three  inductively  tuned  sets 
of  three  series  connected  transducers.  This  multiple- 
inductor  scheme  results  in  both  low  and  uniform  insertion 
loss  among  channels  as  illustrated  in  Figure  9.  This  effi¬ 
cient  scheme  is  especially  necessary  in  the  wide  percentage 
bandwidth  (261),  high  frequency  case.  Here  a  mean  insertion 
loss  of  25.8  dB  with  a  standard  deviation  of  only  0.75  dB 
was  attained.  The  strongest  spurious  tone  for  the  high 
frequency  bank  was  57.5  dB  below  the  desired  tone,  with  only 
about  10%  of  the  spurious  tones  worse  than  63  dB  below  the 
carrier.  Excellent  results  were  also  achieved  for  the  low 
frequency  bank.  The  average  insertion  loss  of  19.8  dB 
(including  1 . 5  dB  due  to  the  switch)  can  be  compared  to  the 
individual  filter  loss  of  14  dB.  A  5  dB  improvement  over 
direct  power  division  is  achieved.  In  addition,  for  the 
low  frequency  bank,  the  worst  case  spurious  was  61.5  dB 
below  the  carrier  with  only  about  3%  of  the  spurious  tones 
stronger  than  63  dB  down. 

/ 

Some  comments  on  these  spurious  levels  are  in  order. 

The  SAW  filters  are  in  fact  capable  of  suppressing  the 
spurious  signals  63-65  dB  as  shown,  for  example,  in  Figure  8. 
In  addition,  under  ideal . conditions ,  the  switches  can  achieve 
70  dB  isolation,  particularly  for  non-adjacent  channels. 

(Note  that  two  switches  were  investigated:  an  approximately 
2  nanosecond  swi-vch  having  an  insertion  loss  of  16.5  dB, 
and  a  250-400  nanosecond  switch  having  an  insertion  loss 
of  1.5  dB)  .  However,'  in  the  finished  filterbanks,  packaging 
and  adjacent  channel  switch  leakage  (possibly  also  due  to 
packaging)  resulted  in  the  levels  quoted. 
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Achieving  high  performance  with  the  SAW  filters  requires 
not  only  proper  design  but  also  careful  fabrication.  Although 
the  transducers  require  lines  and  spaces  as  narrow  as  0.8 
micrometer  together  with  precisd  placement,  they  can  be 
fabricated  with  straightforward  techniques.  Such  fine 
structures  do  require  a  clean  room  environment  and  careful 
work  and  attention  to  cleanliness.  Optical  photoreduction 
from  masters  which  are  10  times  the  final  transducer  dimen¬ 
sions  has  several  advantages,  including  the  relative  low  cost 
of  the  lOX  masters  compared  to  masters  with  0.8  micrometer 
lines  and  spaces,  and  the  very  long  mask  life  due  to  lack  of 
direct  contact.  As  with  any  very  high  resolution  technique, 
monochromatic  light  of  the  shortest  wavelength  compatible 
with  the  optics  and  photoresist  must  be  used.  The  stripping 
method,  in  which  the  area  which  is  to  be  metallized,  is 
exposed,  developed  and  stripped,  before  the  evaporation  of 
the  desired  metal,  is  essential  for  such  high  resolution. 
Using  adjustable  magnification  with  direct  projection  fabri¬ 
cation  reduced  the  mean  deviation  from  design  center  fre¬ 
quencies  to  0.06  MHz  and  0.01  MHz  for  the  high  frequency  and 
low  frequency  filterbanks  respectively. 

It  can  be  concluded  that  SAW  filters  offer  outstanding 
performance  in  compact,  contiguous  filterbank  devices  and 
that  series-parallel  interconnection  is  the  preferred  manner 
of  achieving  the  multiplexing  function. 

Output  Section 

The  output  or  mixer  section  is  the  right-hand  block  of 
Figures  4  and  5  (labeled  MIXER  PACKAGE).  The  high  perform¬ 
ance  in  spurious  rejection  and  low  noise  outlined  in  Figure  6 
requires  that  the  output  section,  like  the  rest  of  the  syn¬ 
thesizer,  be  carefully  designed.  Therefore,  it  is  necessary 
to  use  high  quality  doubled-balanced  mixers,  filters  with 
high  out-of-band  rejection,  and  careful  attention  to  power 
levels.  Fortunately,  these  components  are  not  major  factors 
in  the  synthesizer  cost.  , 

Whenever  two  different  frequencies  are  mixed,  sum  and 
difference  signals  plus  the  input  signals,  multiples  of  the 
input  signals,  and  sums  and  differences  of  one  input  signal 
with  multiples  of  the  other  input,  all  are  present.  The 
original  choice  of  the  frequencies  to  be  generated  and  mixed 
was  made  using  standard  techniques  for  calculating  and  mini¬ 
mizing  the  mixer-generated  spurious  signals.  While  it  is  not 
possible  to  shift  the. frequencies  so  that  all  possible 
spurious  tones  are  outside  the  desired  band,  it  is  feasible 
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to  place  one  p£  the  input  signals  so  that  its  higher  order 
undesired  products  are  all  outside  the  band  of  interest, 
and  set  this  input  approximately  20  dB  stronger  than  the 
other  input.  The  weaker  input  can  then  have  its  higher- 
order  spurious  fall  in-band,  but,  as  they  will  be  extremely 
weak,  they  will  be  adequately  rejected  by  a  balanced  mixer. 
This  was  the  technique  followed  here.  It  is  illustrated 
in  Figure  10,  which  shows  the  most  significant  spurious 
frequencies  at  the  doubler  and  each  of  the  two  mixers.  The 
input (s)  and  the  desired  output  for  each  stage  is  given  by 
the  solid  blocks,  and  the  spurious  signals  by  the  dashed 
ones. 


We  set  up  breadboard  output  section  experiments  in- 
house,  as  well  as  obtaining  a  unit  made  to  our  specifica¬ 
tions  under  AF  contract  by  Communitronics  Ltd.,  Hauppauge, 

NY.  Low  to  moderate-priced  doublers  and  mixers  were  found 
to  have  adequate  performance  provided  that  power  levels  were 
optimized.  As  in  any  case  where  doublers  or  mixers  are  used, 
filtering  is  required.  The  stringent  requirement  to  main¬ 
tain  low  spurious  levels  means  that  the  filters  must  be 
correspondingly  highly  selective.  SAW  filters  would  have 
the  advantage  of  extreme  compactness,  but  at  the  cost  of 
increased  insertion  loss.  A  number  of  low-loss  electromag¬ 
netic  filters  are  available  commercially.  Tubular  filters, 
although  somewhat  bulky  and  expensive  for  a  final  engineered 
version,  have  proved  valuable  for  breadboard  use  in  both  the 
in-house  experiments  and  the  version  built  by  Communitronics. 

The  output  frequencies  obtained  by  mixing  the  1984  MHz 
tone  with  all  allowed  combinations  of  the  two  sets  of  nine 
tones  are  given  in  Table  I.  The  other  frequencies  are 
readily  obtained  by  adding  1  and  2  MHz  to  the  given  values 
corresponding  to  the  use  of  the  1985  and  1986  MHz  tones, 
respectively.  In  the  actual  system  hardware,  logic  circuits 
will  choose  the  proper  combination  of  three  tones  when  a 
given  output  frequency  is  called  up. 

/ 

Summary  and  Conclusions 

The  advantages  of  the  SAW  direct  frequency  synthesizer 
have  been  described.  A  Manufacturing  Technology  contract  to 
reduce  the  size  of  the  present  laboratory  breadboard  to  less 
than  60  cubic  inches  has  been  awarded.  This  contract  will 
draw  upon  the  exverience  of  the  RADC  model  as  well  as  two 
previous  ManTech  contracts  which  improved  the  manufacturabil¬ 
ity  of  SAW  components,  and  develop  a  reasonably  priced  252 
frequency  integrated  JTIDS/TACAN  synthesizer. 
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The  RADC  Deputy  for  Electronic  Technology  held  a  Work¬ 
shop  on  SAW  Frequency  Synthesizers  in  May  1978  to  brief 
industry  on  the  results  of  our  work  and  to  demonstrate  the 
operation  of  the  breadboard  model. 


TABLE  I 


Generation  of  required  frequencies  by  mixing 


high 

freq  bank 

360 

372 

384 

396 

408 

432 

444 

456 

468 

high  freqsi 
doubled 

720 

744 

768 

792 

816 

864 

888 

912 

936 

321 

1585 

1561 

1537 

1513 

1489 

1441 

1417 

1393 

1369 

324 

1588 

1564 

1540 

1516 

1492 

1444 

1420 

1396 

1372 

327 

1591 

1567 

1543 

1519 

1495 

1447 

1423 

1399 

1375 

330 

1594 

1570 

1546 

1522 

1498 

1450 

1426 

1402 

1378 

low 

freq 

333 

1597 

1573 

1549 

1525 

1501 

1453 

1429 

1405 

1381 

bank 

336 

1600 

1576' 

1552 

1528 

1504 

1456 

1432 

1408 

1384 

339 

1603 

1579 

1555 

1531 

1507 

1459 

1435 

1411 

1387 

342 

1606 

1582 

1558 

1534 

1510 

1462 

1438 

1414 

1390 

345 

— 

— 

— 

— ' 

■  — 

1465 

— 

— 

— 

The  above  output  frequencies  are  1984  -  (2  X  high  freq  -  low 

freq) 

Add  one  MHz  for  output  with  1985  MHz  at  the  second  mij^er. 

Add  2  MHz  for  output  with  1986  MHz  at  the  second  mixer. 
Dashed  output  frequencies  are  not  desired. 
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FIGURES 


Figure  1  Single  Channel  Synthesizer 

Figure  2  Photograph  of  a  bank  of  9  SAW  filters  (left) 
and  a  9  X  1  switch  (right)  for  selecting  the 
output  frequency. 

Figure  3  SAW  Filter  Structure 

Figure  4  Frequency  Synthesizer,  Simple  Block  Diagram 

Figure  5  Frequency  Synthesizer,  Detailed  Block  Diagram 

Figure  6  Spurious  and  Noise  Specification  Curves 

Figure  7  Block  Diagram  of  Mode  Locked  SAW  Oscillator 

Figure  8  Insertion  Loss  Versus  Frequency  for  Single 
SAW  Filter 

Figure  9  Insertion  Loss  Versus  Frequency  for  SAW 
Filterbank 

Figure  10  Spurious  Signal  Diagram 


9  TONE  DIRECT  SAW  SYNTHESIZER 


FIG  I 


DIRECT  SAW  MIXER  SYTHESIZER 


FREQUENCY 

FIG  4  HOPPED  OUTPUT 


I  MHz  CLOCK 


SYNTHESIZER  NOISE  AND  SPURIOUS  SPECIFICATION 


FIG  6 


FINAL  VERSION  OF  MLSO 
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Abstract 


A  brief  description  of  Fourier  Spectroscopy  will  be 
presented  as  an  introduction  to  the  description  of  a  novel 
technique  of  potential  application  to  military  surveillance, 
strategic  and  tactical  missions.  The  AFGL  Background 
Optical  Suppression  Scheme  (BOSS)  will  be  described  in 
Hetail  and  its  potential  for  target  detection  in  the  pre¬ 
sence  of  interfering  background  radiation  will  become  evi¬ 
dent.  Experimental  studies  of  the  principle  of  the  technique 
will  be  shown  where  background  suppressions  of  two  orders  of 
magnitude  have  been  obtained.  Modifications  of  the  implemen¬ 
tation  of  the  technique  will  be  briefly  presented.  These 
modifications  allow  background  suppression  to  be  effective 
even  for  structured  background  as  seen  by  an  IR  sensor  system. 
The  technique  allows  spatial,  spectral  and  temporal  dis¬ 
crimination  to  be  accomplished  in  one  sensor  for  enhanced 
mission  performance.  A  simple  implementation  of  the  tech¬ 
nique  will  be  described  for  the  specific  application  of 
auroral  measurements  from  an  aircraft  when  the  measurement 
constraints  are  that  the  sensor  must  look  at  the  target 
through  a  warm  window.  Consequences  of  the  use  of  the  BOSS 
technique  for  this  particular  application  will  be  discussed. 


Ill 


Introduction 


The  essence  of  the  technique  of  Fourier  Spectros¬ 
copy^”^  is  to  take  radiation  of  large  spectral  bandwidth  and 
sinusoidally  modulate  the  intensity  of  each  spectral  compo¬ 
nent  (of  wavelength  •  Xnf  or  wavenumber  a^-l/Xn^  band 

at  a -frequency  proportional  to  and  of  amplitude  propor¬ 

tional  to  its  intensity  I (0^) .  The  radiation  from  all  of 
the  spectral  components  are  simultaneously  detected  by  a 
single  detector  whose  output  is  called  an  interferogram. 

The  instrument  which  is  used  to  produce  this  selective 
modulation  is  an  ordinary  Michelson  interferometer  (or  some 
of  its  modifications)  as  shown  in  Figure  1.  The  radiation 
under  study  is  amplitude  divided  into  two  identical  beams  by 
the  beamsplitter  and  recombined  into  one  beam  after  one  beam 
has  traveled  a  distance  x  greater  than  the  other.  As  the 
optical  path  difference  x  between  the  beams  is  increased  (by 
moving  mirror  Mn )  the  detector  output  yields  the  autocorrela¬ 
tion  function  F(x)  of  the  input  electromagnetic  field.  The 
detector  output  (neglecting  contents  and  numerical  factors) 
is  given  by 


F(x) 


cos  2iTax  da 


(1) 


where  (02-01)  represents  the  radiation  spectral  band  width 
Aa.  The  spectral  power  density  of  the  radiation  1(a)  is  thus 
seen  to  be  given  by  the  inverse  transform 


1(a) 


cos 


2iTax  dx 


(2) 


where  X  represents  the  maximum  path  difference,  produced  by 
the  interferometer,  between  the  two  interfering  beams.  The 
value  X  determines  the  spectral  resolution  6a=l/2X. 

The  interferometer  is  used  for  spectral  measurements 
because  of  its  manifold  advantages  over  conventional  grating 
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or  prism  techniques  .  However,  it  also  turns  out  that  it 
lends  itself  easily  to  a  technique  of  background  suppression^ 
which  we  shall  now  describe. 

Background  Optical  Suppressioh  Scheme 

Let  us  consider  the  two  interferometers  shown  in 
Figure  2  with  their  corresponding  interferograms  shown  on 
their  right.  We  notice  immediately  that  the  interferograms 
are  complementary,  while  the  only  change  as  regards  the  inter¬ 
ferometers,  is  that  the  input  radiation  enters  the  inter- 
feroiheter  through  the  upper  face  of  the  beamsplitter  for  the 
upper  case,  but  the  lower  face  of  the  beamsplitter  for  the 
lower  case.  It  should  be  evident,  that,  if  the  radiation  is 
allowed  to  enter  a  single  interferometer  simultaneously 
through  the  upper  and  lower  faces  of  the  beamsplitter,  (which 
would  correspond  to  adjacent  fie Ids-of- view  of  the  interfero¬ 
meter)  the  two  interferograms  at  the  right  in  Figure  2  would 
be  simultaneously  generated  and,  upon  falling  onto  one 
detector,  would  yield  an  interferogram  which  has  no  struc¬ 
ture.  This  is  tantamount  to  saying  that  no  spectral  informa¬ 
tion  about  the  radiation  is  contained  in  the  output  inter¬ 
ferogram  function. 

Let  us  assume  that  the  input  radiation  is  sky  back¬ 
ground  radiation  which  is  reasonably"*"  uniform  and  enters  the 
interferometer  through  the  upper  and  lower  faces  simultan¬ 
eously.  Because  of  the  complementary  characteristics  of 
their  interferograms,  the  resultant  detector  output  shows  no 
modulation.  Imagine  now  that  a  target  appears  in  one  field- 
of-view  but  not  the  other;  the  signal  generated  by  the  back¬ 
ground  is  a  constant  but  now  the  interferogram  has  structure 
which  is  due  to  the  target  only.  A  Fourier  transformation 
of  this  interferogram  yields  the  spectral  signature  of  the 
target  with  no  spectral  content  due  to  the  background. 

Demonstration  and  Application 

Figure  3  shows  experimental  results’*"'’  obtained  under 
conditions  of  uniform  background  and  adjacent  f ields-of-view. 


*The  interested  reader  can  refer  to  the  references  for  a 
more  rigorous  and  complete  description  of  Fourier  Spectro- 
s  copy . 

+This  suppression  technique  can  also  be  implemented  to  sup¬ 
press  structured  background,  but  this  will  not  be  discussed 
here  due  to  lack  of  space;  see  reference  5* 

++Work  accomplished  under  contract  with  Visidyne,  Inc.,  and 
funded  by  the  AFGL  Laboratory  Director's  Fund. 
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The  upper  left  trace  shows  the  background  radiation  coming 
from  one  field-of-view. only,  while  the  upper  right  trace  is 
from  the  other  f ield-of-view  only.  The  bottom  trace  is  the 
resultant  interferogram  when  the  background  enters  the  two 
f ields-of-view  simultaneously.  The  background  suppression  . 
was  measured  to  be  about  two  orders  of  magnitude;  which 
implies  that  the  target  could  be  fainter  than  the  background 
by  about  100  and  still  geners,te  detectable  structure  in 
the  interferogram. 

We  have  mentioned  above  that  structured  backgrounds  can 
be  treated  as  well  as  uniform  backgrounds  (although  as  must 
be  noted  with  increased  complexity) .  Here  we  propose  using 
the  technique  for  the  case  when  the  interferometer  must  view 
a  faint  or  distant  target  through  a  warm  window,  where  the 
emission  from  the  window  is  significantly  greater  than  that 
from  the  target.®  The  obvious  approach  to  suppressing  the 
window  spectral  emission  is  to  operate  the  interferometer  in 
the  dual-beam  mode  where  the  window  emission  enters  simul¬ 
taneously  the  two  f ields-of-view  of  the  interferometer, 
while  the  target  is  seen  through  only  one  f ield-of-view. 

One  example  involving  the  use  of  a  dual  beam  interfero¬ 
meter  is  the  spectral  measurement  of  infrared  auroral  emis¬ 
sion  from  an  aircraft  as  illustrated  in  Figure  4.  The  air¬ 
craft  can  be  flown  at  sufficiently  high  altitude  to  avoid 
much  of  the  absorption  and  emission  of  the  lower  atmosphere, 
but  even  at  high  altitudes  the  emission  of  the  atmosphere 
in  the  infrared  is  significant.  A  further  complication  is 
added  by  the  emission  of  the  window  through  which  the  inter¬ 
ferometer  views  the  aurora.  The  interferometer  operating 
dual  beam  only  measures  the  difference  between  what  is  view¬ 
ed  by  the  two  beams.  Since  both  beams  view  the  window  and 
the  atmosphere,  the  resultant  interferogram  contains  no  mod¬ 
ulation.  When  an  auroral  arc  enters  one  beam  there  is  a 
difference  and  the  interferometer  signal  is  an  interferogram 
of  only  the  auroral  emissions. 

The  measurement  of  auroral  emission  is  just  one  ex¬ 
ample  of  the  potential  uses  of  the  dual  beam  background 
suppression  interferometer.  In  geneial,  the  technique  can 
be  used  anytime  a  target  can  be  spatially  isolated  against 
a  uniform  background  and  could  certainly  be  used  to  measure 
aircraft  and  rocket  emissions  against  warm  backgrounds. 
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FIGURE  1 


EXPERIMENTAL  RESULTS  ILLUSTRATING  THE  EFFECTIVENESS  OF  THE  AFGL 
BACKGROUND  SUPPRESSION  SCHEME 
(Photograph  Courtesy  of  Visidyne) 
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FIGURE  3 


118 

I 


SPECTRAL  MEASUREMENT  OF  INFRARED  AURORAL  EMISSION  USING  A 
DUAL-BEAM  BACKGROUND  SUPPRESSION  INTERFEROMETER 
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ABSTRACT 

An  adaptive  control  methodology  is  developed  based  on 
optimal  principles  using  spectrum  estimation  procedures  of 
complex,  scatterer  resolved,  lohg  irange  target  returns  for 
purposes  of  precision  pointing  control  and  selected  target 
scatterer  illumination.  Aim  point  selection  of  scatterers 
would  assist  in  target  identification  and  discrimination 
operations.  It  is  the  intent  of  the  adaptive  controller/ 
spectrum  estimator,  developed  here,  to  assist  in  such 
operations  in  such  a  manner  that  range  transit  time  does  not 
‘  ultimately  limit  adaptive  pointing  control  and  scatterer 
selection  on  target,  and  that  the  scattering  spectrum  inherent 
in  complex  structured  targets  are  fully  used  as  a  new  dimen¬ 
sion  of  information  for  pointing  system  control  purposes. 

The  approach  incorporates  autoregressive  spectrum  estimations 
algorithm  techniques  for  scatterer  induced  spectrum  resolu¬ 
tion  and  coherent  spectrum  detection.  Adaptive  principles 
are  based  on  scatterer  spectrum  selection  criteria  with 
spectrum  width  convergence  serving  as  target  loop  closure 
logic.  -The  control  algorithm  is  modeled  in  an  optimal  sense 
through  the  selection  of  spectrum  and  motion  observables 
which  demonstrate  enhanced  sensitivity  to  scatterer  conver¬ 
gence  and  spectrum  noise  discrimination.  An  adaptive  optics. 
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complex  structure  target  example  is  considered  with  energy 
convergence  as  a  performance  figure-of-merit  at  a  selected 
scatter.  This  example  will  use  experimentally  obtained 
scatterer  data,  and  an  asstamed  target  structure  and  active 
pointing  system  parameter  set  which  is  representative  of 
current  long  range  identification/discrimination  mission 
operations  of  interest  to  the  Air  Force. 


OVERVIEW  OF  APPROACH  AND  DEVELOPMENT 

This  analysis  addresses  active  illuminating  systems  and 
complex  targets.  That  is,  targets  characterized  as  having 
structure  larger  than  the  beam  size  in  the  target  space  and 
scattering  surface  texture  which  is  varied  in  nature  locally 
across  the  target,  with  variations  which  include  specular 
and  diffuse  textures  relative  to  the  illuminating  wavelength 
and  correlation  length.  The  illuminating  system  is  assumed 
to  be  adaptive  both  in  transmitter  wavefront  phasing  and  in 
pointing.  That  is,  wavefront  focusing  or  phasing  and  energy 
•centroid  steering  or  positioning  are  assumed  to  be  available 
in  real  time  for  control  in  the  target  space.  Control 
criteria  is  energy  maximization  at  a  resolvable  target 
scatterer,  with  scatterer  being  both  identifiable  and  select¬ 
able  by  the  system  controller.  The  study  results  will,  show 
that  backscatter  field  spectrum  narrowing  provides  a  good 

convergent  criteria  which  inherently  leads  to  energy  maxiraiza- 

/ 

tion,  and  that  spectrum  correlation  can  aid  in  scatterer 
selection  for  aimpoint  control. 

Because  complex  target  returns  behave  randomly  at  the 
receiver  (random  in  regards  to  their  point  of  origin  and 
intensity  due  to  constructive  interaction  between  specular 
and  diffuse  returns)  an  entirely  r^dom  active  system  is 
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considered ,  that  is  a  random  pulse  transmitter  is  part 
of  the  modeled  active  system  with  spatial  coherence  and 
adaptation  being  the  system  control  parameters .  The 
adaptive  controller  is  based' on  spectrtun  estimation  via 
autoregressive,  methods ,  with  spectrum  frequency  narrow¬ 
ing  as  the  prime  loop  closure  logic  for  energy  convergence  on 
a  selected  target  specular  locality,  or  on  a  finite  extent 
target  diffuse  region.  A  second  motivation  for  selecting  a 
random  signal  system  approach  is  the  freedom  from  range  and 
frequency  ambiguity  due  to  the  random  waveform.  This 
inherently  will  allow  the  adaptively  controlled  system  to 
achieve  maximum  target  scatterer  resolution  and  energy 
maximization  without  waveform  ambiguity  limitations.  This 
feature  optimizes  the  employment  of  the  super  resolution 
properties  of  autoregressive  spectrum  estimation  techniques 

9 

♦  . 

We  will  first  develop  the  received  signal  scattering 
model,  the  system  correlation  processor,  and  the  controller 
architecture  applicable  to  adaptively  phase  driven  random 
pulsed  coherent  systems.  This  will  be  followed  by  a  perfor¬ 
mance  assessment  of  super  resolution  techniques  simultaneously 
being  employed  as  adaptive  controller  and  signal  parameter 
estimation. 
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2.  MODELING  -  SIGNAL,  SYSTEM,  CONTROLLER,  ADAPTIVE  CRITERIA 


Since  the  goal  of  this  development  is  high  resolution 
of  target  structure  and  scattering  centers  as  well  as 
adaptive  illumination  (with  maximum  energy  convergence)  of 
a  selected  scattering  center,  the  illiiminating  system  chosen 
is  a  pseudo  random  waveform  pulsed  radar  with  w'avefront  phase 
control  implemented  on  a  pulse-to-pulse  basis.  The  random 
radar  is  selected  due  to  its  freedom  from  range  and  velocity 
ambiguities ^ ,  permitting  maximization  of  target  resolution  ‘ 
both  in  range  (slant  range)  and  velocity  (cross  range). 
Pulse-to-pulse  wavefront  phasing  is  selected  for  energy 
convergence  to  a  single  target  region  or  scattering  center' 
for  enhanced  target  resolution  in  support  of  high  resolution 
identification  missions. 


A  received  pulse  x(t')  at  time  t’  can  be  related  to  a 

family  of  proceeding  pulses  x(t'  -  kT) ,  k  =1,  L  (defined 

as  the  pulse  sequence  5x(t')  7  ),  all  having  a  random 

i-  3  k ,L 

waveform  but  a  specified  wavefront  phase  distribution,  by  an 
Lth  order  autoregressive  process  given  by  the  follovring 
difference  equation. 


x(tV) 


a 


(L) 

X  (t 


-  kT)  +  n(t') 


(1) 


k=l 
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where  T  iis  the  pulse  period,  n(t')  is  a  stationary  Gaussian 

•  ■  'i  .  ’ 

white  noise  process  with  zero  mean  and  variance 

associated  with  the  signal  sample  x(t').  The  digitized 

version  of  this  form,  assuming  sample  time  spacing  T, 

sequence  of  available  time  samples  and  current  time 

L 

signal  sample  Xj^,  becomes 


(L) 

^i-k 


which  defines  the  current  signal  sample  in  terms  of  the  L 
formet  samples  (L  lags  for  the  autoregressive  process). 
Coefficients  will  be  defined  later  in  terms  of  the 

receiver  structure  as  a  sampled  sequence  correlator. 

That  is,  for  enhanced  resolution  and,  sensitivity,  the  receiver 

•  • 

will  be  modeled  as  a  pulse  correlator  which  references  a  sample 
of  x(t)  that  has  been  stored  by  the  system  at  the  time  of. 
transmission  for  purposes  of  coherent  detection.  Returning 
to  the  time  analog  format  (Eq  (1)),  each  pulse  experiences 
a  transit  time  lag  Tj^(t)  =  "  ^k^'  where  is  the  lag  rate 

at  t,  the  time  referenced  to  pulse  transmission  which  is 
related  to  the  pulse  arrival  time  t'  as. 
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t*  ="t  -  Tk(t)  =  (1  +  C^)  t  -  T,^  (3) 


(4) 

The 


(5) 

r(t)  =  ^x(t  -  Tj(t))}  =  (1  -  Cj)  t  -  ' 

.  L  L 

"  A'  CORKELATOR-COHEPENT  DETECTION; 

The  correlator  output  can  be  modeled  as  a  low  pass  ‘ 
filtered  output  (the  expected  value  of  the  product  x(t)  r  (t) ) 
having  the  fOna, 

E  (x(t)  r(t)l  =  R  _(t) 

» 

.  129. 


giving , 


(t)  =  ^  X  ((1  +  Cj^)t  -  Tk  -  kT)  +  nCt) 

k 


noting  the  initial  incremental  range  of  k  =  1,  L. 
correlator  will  reference  some  jth  pulse  in  the  stored 
sequence  {x(t)j^  yielding  the  following  reference 
sequence 


[stored  sample  of  x(t)]  -►  r(t),  where 


L,L 

=  Hi  -  (VTj)i 

k.j  ■  .  ■ 

+  E  In(t)  r(t)]  (6) 


where  is  the  autocorrelation  of  x(t).  In  the  sampled 
system  where  both  signal  and  reference  have  been  digitized 

i 

as  a  binary  data  set  with  sample  intervals  of  length  T  and 
where  reference  signal  correlation  with  transmitted  wave¬ 
form  noise  (n(t))  gives  a  non  zero  value  only  for  the  data 
sample  occurring  at  time  t,  the  correlator  output  can  be 
modeled  as , 

.  ■ ,  ■  ^  ■  ;  ■ '  ■  '  .  L  ^ ' 

E  (L) 

i  r  : 

•  * 

+  E  [njL  r^l 


-  (7) 


B.  RECEIVED  NOISE  ENERGY-MAXIMIZATION  CRITERIA 

■  If  we  return  to  Eq  (2),  the  digitized  L  lag  format 
for  current  time  sample  x^  (at  t  -  t'  =  i  T) ,  we  can  define 
ao  =  “1/  giving  for  the  received  noise  energy  at  current 
time  increment  i. 
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.5^—1 

Hi  =  2_.  aj, 

k=l 


(L)  Y-' 

c  xi-k  -  =  L 


(L) 

ajc  Xi_k 


k=o 


Thus  noise  energy  conservation  demonstrated  by  the  random 

waveform  radar  can  be  defined  by  the  correlator  stored 

reference  sample  set  -T plus  the  current  time  sample 

L  . 

in  terms  of  total  noise  at  the  current  time  t  by  Eq,  (7). 

A  condition  of  phasing  the  signal  pulses  at  the  target  which 
maximizes  the  pulse  summing  in  Eq,  (8)  can  readily  serve  as 
a  target,  illumination  control  criterion  for  energy  maximization 
at  a  single  target  scattering  point.  Thus  phasing  adaptation 
is  based  on  ■ 

L  ■ 


nil 

maximization 


(L) 

2_ _ I  ^k  Xi_k 

k=o 


✓ 


(9) 


for  a  random  waveform  radar  using  coherent  pulse  detection 
based  on  the  sequence  I  3C  1  , 
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•> 

C.  COMPUTATION  RELATIONS  AND  NOISE  VARIANCE; 

Since  correlates  only  with  (current  time 
sample),  multiplying  Eq.  (8)  by  n^  and  taking  the  expected 
value  gives. 


E  [  n^^]  =  E  [ni  Xi  )  = 


(10) 


This  allows  us  to  define,  in  a  mean  squared  error  sense,  from 
Eq.  (7), 

■  ,  L 

=  E  [  x(t)  r(t)  ]  -  ®i^'  ^xx 

i=l  * 


=>xx(o)  ^ 

i=l 


1=0 


(11) 


where  the  coefficients  (noting  =  ‘^1)  can  now  be 

L 

defined  in  terms  of  the  sample  correlation 


N-i 


k=l 


(12) 
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and  the  following  recursion  relationships,  bearing  in  mind 
that  aj^  denotes  the  ith  autoregressive  coefficient  of  an  Lth 
order  autoregressive  process,  i  <  L, 


RECURSION  RELATIONSHIP: 


a^^^V=  a^  ;  i  =  1,  . . . ,  L-1  (13) 


(L) 


with  the  final  coefficient  a^^  •  satisfying 


FINAL  COEFFICIENT: 


L“1 


R(L) 


(L) 


-E-‘ 


(L-1) 


R  (L-i) 


i=l 


i=l 


Ii(i) 


(14) 


*  This  sequence  is  worked  out  in  Section  3  of  this  paper  where 

a  performance  assessment  of  the  spectrm  estimator/controller 

2 

is.  made.  Now  that  we  have  the  noise  variance  and  control 

criterion  n^]  (Eqs.  (11)  and  (9),  respectively)  modeled  in 
max 

terms  of  the  sequence  correlator,  we  can  define  the  target 
scattering  spectrum  in  terms  of  and  a  scattering 
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correlation  time  b“^,  relatable  to  the  coherent  integration 
interval  NT,  or  the  number  of  total  scattered  pulse  sequences 
which  can  be  coherently  summed  ^ ^ ,  b""^  is  driven  by  the 

doppler  spread  in  the  backscattered  return,  being  different 
for  glint,  diffuse  and  self  convolved  adaptive  control  spectral 
scattering  regions.  An  examination  of . the  spectrum  profile 
in  Figure  1  should  clarify  this  point.  Glint  (see  Note  1)’  or 
specular  returns  generally  originate  from  discontinuous  and/or 
smooth  (in  terms  of  the  transmitting  wavelength)  target  regions, 
where  the  latter  case  involves  target  geometry  in  terms  of  the 
directionality  (towards  to  receiver  line-of-sight)  of  the 
local  surface  normal.  These  returns  generally  remain  coherent 
with  system  stored  transmitter  references,  demonstrating 
large  bg“^  magnitudes  (or  narrow  spectral  peaks.  Thus  the 
glint  scattering  spectrum  can  be  modeled as, 

Sgg  =  ^  bg  o„2/  (w2  +  bg2)  US) 


Note  1.  The  definition  of  specular  or  glint  as  used  here  is 
based  on  the  optical  or  near-optical  region  of  the  spectrum  where 
returns  from  wavelength  smooth  surfaces  preserve  illuminating 
spatial  coherence. 
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The  diffuse  (surface  texture  regions  rough  or  on  a 
comparable  scale  to  the  illiiminating  wavelength)  generally 
represent  extended  target  regions  of  the  dimensions  of  the 
illuminating  be^  size  in  the  target  space.  The  backscatter 
correlation  time  b^j”^  is  generally  quite  short,  giving  a 
much  broader,  reduced  magnitude  spectrum  for  the  diffuse 
returns,  with  typically 

< 

b^  >  10  bg  ;  thus, 

|Sdd(o)UlSgg(o)  1/10  (16) 

The  self  convolved  generated  spectrum  eminates  from 
spectrum  modulation  due  to  adaptively  driven  illxamination  of 
selected  target  localized  regions  convolved  with  the  target 
illviminated  silhouette  which  is  represented  by  the  sampled 
signal  sequence,  that  is,  the  time  history  of  illxamination 
coivered  by  the  system  for  the  L  history  samples.  The  control 
algorithm  which  was  found  adequate  was  based  ph  a  simple 
monotonic  stepping  routine,  where  wavefront  convergence  or 
phase  tapering  implemented  a  series  of  corrections  in  one 
directional  sense,  then  returned  to  a  nominal  profile  (contained 
within  or  at  the  start  of  this  sequence)  in  anticipation  of 


adaptive  performance  assessment.  This  assessment  was  based 
on  Eq.  (9)  criteria,  that  is  maximization  of  the  received 
random  (noise)  transmitter  energy,  for  a  ttial  sequence  of 
wavefront  profile  changes.  Performance  assessment  involved 
the  simple  test  of  n^  increase  trend  sensing  over  a  L-N 
trial  sequence,  where  N  approached  L-1  when  energy  increase 
trends  were  confirmed.  Wavefront  profile  stepping  rates 
were  fixed,  being  nominally  less  than  (10  T)"^  for  stable 
cdnttol  o^nisnce  1 

.  i 

Thus  f control 

The  self . convolved  spectrum  was  adequately  modeled  as 

See  ■=  (w2  -  w^2  )  -  bj,2w2],  '  ■  .  ■  (18) 

w,2  “  2ir  f control  , 

and  b_Aj  b_,  since  stepping  control  preser^res  the  spatial 
coherence  with  the  transmitted  pulse  train,  and  convergence 
generally  hunts  target  regions  which  are  specular  in  nature, 
that  is  preserves  the  spatial  phase  of  the  return.  The 
above  spectral  modeling  is  analogous  to  a  Markov  first  order 
random  process  for  the  specular  (Sgg)  and  diffuse  (S^j^) 

136 


spectral  regions,  and  a  Markov  second  order  process  for 

« 

the  self  convolved  region.  ‘ 

The  estimator/controller  approach  modeled  above, 
namely  (a)  Spectrum  Estimator,  Egs.  (11) ,  (13)  through  (16) , 
and  (18) ;  (b)  Controller  and  Criteria,  Eqs.  (8),  (9)  and 
(17) »  are  tested  in  a  system  example  ih  Section  4. 

3.  PERFORMANCE  ASSESSMENT  OF  SPECTRUM  ESTIMATOR/ADAPTIVE 

CONTROLLER  APPROACH 

We  shall  assess  the  performance  of  the  autoregressive 
spectrum  estimator/adaptive  controller  modeling  approach  by 
considering  a  high  resolution  system  example.  We  select  a 
laser  radar  system  candidate  for  this  demonstration  since 
such  systems  inherently  support  high  resolution  capabilities 
but  also  suffer  from  good  control  procedures  and  criteria 
due  to  target  backscatter  interaction  with  beam  convergence 
controllers. 

'  .  The  candidate  system  is  a  random  waveform  pulsed  radar, 
nominal  pulse  repetition  and  digital  sampling  rates  of  1  MHz, 
operating  wavelength  in  the  near  infrared  at  10.6  micrometers. 
Adaptive  control  involves  phasing  the  outgoing  wavefront  of 
a  sequence  of  pulses  to  maximize  illumination  at  a  selected 
target  scattering  point  and  to  converge  the  beam  size  in  the 
target  space  such  that  a  reduced  target  region  is  illuminated. 
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the  reduction  primarily  limited  by  the  corrected  diffraction 

'  .  i  .  ' 

perfoirmance  of  the  system  antenna  (Cassegrain  beam  expander) 
plus  propagation  path  anomalies  (atmospheric  turbulence  and 
beam  disturbance  inducing  factors) .  Corrected  near-diffraction 
performance  is  experience  htere  because  only  focus  control 
(outgoing  wavefront  phase  profile  on  a  pulse-to-pulse  basis) 
is  implemented.  High  order  wavefront  correction  would 
demonstrate  close  to  vacuum  diffraction  beam  convergence  in 
the  target  space.  The  system  parameters  are  given  in  Table  1. 
The  autoregressive  computing  relations  are  listed  in  Table  2 
for  the  first  3  orders.  System  convergence  is  displayed  in 
Figure  2.  The  convergence  time  profile  is  compared  with  an 
earlier  reported  estimation  technique ^ .  The  same  intensity 
maximum  is  reached  at  the  target  normal  glint  point  for  both 
the  autoregressive  approach  and  the  maximization  estimation 
approach  (diescribed  in  Reference  1)  but  the  fomer  is  at 
least  faster  by  a  factor  of  100.  Indicated  convergence  time 
is  less  than  1  millisecond  for  the  autoregressive  method. 

This  is  primarily  attributable  to  the  non  periodic  nature 
of  the  pulsed  radar  system  configured  here  as  well  as  the 
random  phase  stepping  method.  Random  stepping  is  not  limited 
to  any  system  carrier  frequency  and  can  exhibit  rapid  control 
convergence,  providing  the  convergence  trends  are  estab¬ 
lished  early  in  loop  closure  execution. 
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CONCLUSION 


Spectrum  estimation  and  system  convergence  can  be 
achieved  based  on  pulse  sequerlcc  energy  maximization  logic 
and  autoregressive . estimation  methods  for  non  periodic 
control  procedure  systems.  The  |iulsed  radar  example  con¬ 
sidered  here  has  verified  this  point.  Random  stepping  and 
randbm  signal  radar  systems  using  these  principles  potentially 
can  demonstrate  enhanced  target  irradiance  and  potential 
resolution  increases  over  non  adaptive  or  periodic  controlled 
system  concepts. 
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TABLE  1. 


SYSTEM  PARAMETERS 


^vffc.'<v^-ertp/r 
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TABLE  1. 


SYSTEM  PARAMETERS  (cont) 
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TABLE  2.  COMPUTATION  RELATIONSHIPS  ^ 
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TABLE  2. 


COMPUTATION  RELATIONSHIPS  (cont) 
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FIGURE  1.  NOMINAL  TARGET  SCATTERING  SPECTRUM 

Note  contributions  froni  specular  (glint) 
diffuse  and  adaptive  seif  convolved 
scattering  phenomena' 
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POSITION  OF  FOCUSING  OPTICAL  ELEMENT 
IN  METERS 


FIGURE  2.  FOCUS  CONTROL  CONVERGENCE 


Note  the  comparison  of  the  autoregressive 
approach  (the  stepiped  curves)  versus  an 
earlier  modeling  of  a  estima tor/controller 
based  on  energy  maximization  techniques 
(see  Reference  1) 
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SPECTRUM  POWER,  RELATIVE  MAGNITUDE 
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FIGURES.  SPECTRUM  CONVERGENCE 


to 


Note  spectrum  convergence  as  beam 
phases  to  single  target  scatterer. 

Selected  scatter  was  geometrically 
normal  for  this  example. 
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spatial  and  Temporal  Coding  of  GaAs  Lasers  for 
a  Laser  Line  Scan  Sensor 

Abstract 


Line  scanning  laser  sensors  which  measure  range  and  reflectance  are 
currently  being  developed  to  be  used  as  the  front  end  of  real-time  image 
processors.  This  paper  establishes  the  space-time  combination  of  intensity 
modulation  and  scan  pattern  for  N  GaAs  laser  sources  which  results  in 
optimum  system  performance.  A  noise  model  based  upon  the  Poisson  point 
process  behavior  of  a  direct  (optical  intensity)  detector  is  developed. 
Estimation  theory,  based  upon  the  Poisson  statistics  of  the  detector, 
is  used  to  develop  the  maximum-likelihood  processor  for  both  range  and 
reflectance.  Performance  expressions  are  presented  to  relate  the 
variance  of  the  range  and  the  reflectance  estimates  to  the  line  scan 
system  parameters;  and  the  variance  of  both  is  shown  to  be  proportional 
to  the  detected  signal  energy.  Various  combinations  of  spatial  and 
temporal  coding  of  the  N  laser  sources  to  separate  the  return  signal, 
energy  from  each  ground  resolution  cell  while  increasing  the  returned 
energy  per  cell  are  evaluated.  Based  upon  the  variance  of  the  range 
and  reflectance  estimates,  the  optimum  spatial-temporal  coding  com¬ 
bination  is  chosen.  It  is  shown  that  if  a  large  maximum  unambiguous 
range  is  not  required,  the  best  performance  possible  is  achieved  with 
N  sources  scanned  in  parallel,  each  source  being  sinusoidally  modulated. 

If  a  large  max;imum  unambiguous  range  is  required,  however,  the  best 
performance  is  achieved  with  N  sources  scanned  in  parallel,  each  source 
using  sinusoidal  modulation  which  is  On-Off-Keyed  by  a  pseudo-noise 
(PN)  code. 


Introduction 


Line  scanning  laser  sensors  are  currently  being  developed  to  be  used 
as  the  front  end  of  real-time  image  processors.  Properly  coded  intensity 
modulation  of  the  active  illumination  (scanning  lasfer  beam)  provides  the 
capability  to  record  both  images  (reflectance  data)  and  three  dimensional 
characteristics  (slant  range  data) .  Range  gating  might  be  used  in  the 
classical  radar  sense  for  clutter  rejection  in  order  to  increase  the 
signal  return  from  a  particular  target  in  the  scene  or  to  produce  height 
profiles  (croSs-sectlonal  boundaries).  The  target  height  and  reflectance 
data  could  then  be  processed  by  pattern  recognition  algorithms  for  tar¬ 
get  cueing  and  threat  assessment. 

Future  operational  constraints  require  small,  light,  efficient,  and 
rugged  laser  line  scanners.  A  very  favorable  source  candidate  is  the 
s^lconductor  GaAs  laser  diode  which  is  very  small  and  lightweight.  Un¬ 
fortunately,  the  GaAs  laser  is  also  a  low  power  device,  thus  requiring 
the  use  of  many  laser  diodes  to  produce  the  same  performance  achievable 
with  a  larger,  more  powerful  solid  state  or  gas  laser. 

The  problem  is  to  determine  an  efficient  way  to  combine  N  (an  Integer 
number)  GaAs  laser  beams  for  estimating  range  r  and  reflectance  p.  Un¬ 
fortunately  it  is  impractical  to  combine  N  GaAs  laser  beams  into  one 
beam;  the  optical  problems  involved  are  many  and  complex,  and  large 
single  substrate  arrays  of  GaAs  diode  sources  are  not  yet  within  the 
current  state-of-the-art. 

The  use  of  N  separate  beams  appears  to  be  the  best  approach.  How¬ 
ever,  in  order  to  keep  from  blurring  the  image  and  from  smoothing  the 
three-dimensional  surface,  it  is  necessary  to  separate  the  reflected 
energy  from  each  ground  resolution  cell,  or  pixel.  Hence,  coding  the 
modulation  of  the  N  laser  beams  is  necessary. 


In  order  to  gain  the  full  advantage  of  the  multiple  source  system, 
the  N  transmitter-receiver  pairs  (or  channels)  must  not  interfere  with 
each  other.  Interpixel  or  Interchannel  Isolation  may  be  achieved  by 
making  the  transmitter-receiver  pair  fie Ids -of -view  (FOV)  spatially 
orthogonal,  i.e.,  the  pixel  subtended  by  the  FOV  of  transmitter- 
receiver  pair  'm'  cannot  overlap  the  FOV  of  transmitter-receiver 
pair  ’j'.  Interchannel  crosstalk  may  also  be  minimized  by  constraining 
adjacent  source  modulating  signals,  Sjjj(t)  and  s.(t),  to  be  temporally 
orthogonal,  i.e.,  have  zero  cross-correlation.  It  will  be  shown  that 
the  optimum  spatial-temporal  coding  combination  utilizes  both  spatially 
and  temporally  orthogonal  signals . 


In  the  sections  to  follow,  a  brief  summary  of  the  statistical 
model  for  optical  detectors  which  operate  in  the  quantum  limited  re¬ 
gime  is  presented.  This  model  is  then  used  to  obtain  near-optimum 
estimators  (post  detection  processors)  of  both  pixel  reflectivity  and 
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slant  range  from  the  aircraft.  The  Cramer-Rao  bound  is  included  as  an 
indication  of  estimation  performance  for  the  various' modulation  candi¬ 
dates.  A  number  of  laser  modulation  formats  are  presented  and  their 
relative  performance  compared. 

System  Model 

This  section  developes  the  mathematical  tools  necessary  to  compare 
the  performance  of  the  various  laser  modulation  formats.  To  begin  with, 
a  noise  model  based  upon  the  Poisson  point  process  behavior  of  a  direct 
(optical  intensity)  detector  is  developed.  Estimation  theory,  based 
upon  the  Poisson  statistics  of  the  detector,  is  then  used  to  develop  the 
maximum-likelihood  processor  for  both  range  and  reflectance.  In  the 
final  part  of  this  section,  performance  expressions  are  developed  to 
relate  the  variance  of  the  range  and  the  reflectance  estimates,  and 
other  performance  criterion,  to  the  laser  line  scan  system  parameters. 

Optical  Detector  Statistical  Model.  Due  to  the  uncertainty  in  photon 
arrival  times,  the  quantum  nature  of  the  photoelectron  generation  process, 
and  the  random  nature  of  background  noise  and  dark  current,  the  output  of 
an  optical  detector  is  a  stochastic  process.  Conditioned  on  knowing  the 
signal  component  of  the  complex  envelope  of  the  incident  optical  field, 
the  detector  output  is  conditionally  an  inhomogeneous  filtered  PoissOn 
point  process  with  the  rate  function  [1], 

>^(t)  =  ^Sa  !Ug(?.t)l2d?  + 

Pj[  E[IUb(f,t)I2]df  +  Xd  (1) 


where: 

average  number  of  electrons  generated  per  unit  of 
incident  optical  energy 
detector  area 

complex  envelope  of  the  incident  electromagnetic 
field  due  to  a  known  optical  signal 
complex  envelope  of  the  electromagnetic  background 
noise  field 

dark  current  rate  function 

expected  value  (ensemble  average) 
operator 


3  = 


A  « 

U  = 
s 


E  = 
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It  is  well  known  that  the  probability  of  detecting  m  photons  in  the 
time  interval  [0,T],  that  is  N(T)  =  m,  is  given  by 

P[N(T)  =  m]  =  e  ;  m  =  0,1,2,  ... 

m .  .  . 

where 

y  =  r  X(a)da 
Jo  ¬ 
in  addition,  conditioned  on  knowing  U  ,  the  joint  probability  density  of 
the  ordered  event  times  {t^}  and  the  number  of  events  N(T)  in  an  interval 
[ 0 , T ]  is  given  by , 

=  X(t^)  '  (2) 


This  is  a  complete  statistical  characterization  of  the  photo- 
detector  output.  The  joint  probability  density  function  (pdf)  given  in 
Eq  (2)  will  be  the  basis  for  the  parameter  estimation  which  follows. 

Parameter  Estimation.  The  principle  behind  the  measurement  of 
range  and  reflectance  by  the  laser  line  scan  system  is  very  similar  to 
that  of  a  classical  radar.  As  with  a  radar,  the  time  delay  (or  phase 
difference)  between  a  transmitted  signal  and  its  return  is  proportional 
to  the  target  range.  In  addition,  the  return  signal  amplitude  is  pro¬ 
portional  to  the  target  reflectance. 

All  measurements  of  the  actual  values  of  the  delay  t  and  the  re¬ 
flectance  p  will  be  degraded  by  system  noise,  and  thus  will  only  be 
estimates  of  the  true  values.  These  estimates,  denoted  t  and  p,  are 
based  upon  the  measurements  made.  In  the  case  of  well  designed  optical 
detector,  the  observations  are  the  ordered  event  times  (t^}  of  detector 
output  and  the  nimber  of  events  N(T)  in  an  interval  [0,T].  The  desire 
is  to  use  statistical  tools  available  to  derive  a  signal  processor  such 
that  E[t]  =  T  (the  estimate  has  zero  bias)  and  var  [t-t]  (a  measure  of 
how  close  the  estimate  t  is  to  the  real  value  t)  is  small,  with  similar 
expressions  for  p  and  p.  Since  no  a  priori  information  about  the  random 
parameters  t  and  p  is  known,  maximum  likelihood  (ML)  estimation  is  used. 
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Given  the  observation  Z,  the  ML  estimate  of  the  random  parameter  A, 
is  defined  by 


Maximize  An  f-i. 
by  Choice 
of  A  ^ 


(z|a) 


[2,p.6]  where  f^j^CzjA)  is  the  conditional  pdf  of  the  observation  Z, 

conditioned  on  knowing  the  parameter  A.  For  convenience,  the  subscripts 
will  be  dropped  from  the  traditional  notation  for  pdf's  and  N  will  be 
used  to  represent  N(T)  =  m.  From  Eq  (3)  it  follows  that  a  sufficient, 
though  not  necessary,  condition  for  the  ML  estimate  can  be  found  from 


3  iln[f(Z[A)l 
3A  A 


If  the  signal  dependent  portion  of  Eq  (1)  is  replaced  by  ps(t) 

and  the  remaining  portion  is  replaced  by  X  ,  then  Eq  (1)  becomes 

n 


X(t)  =  ps(t)  +  X^  (5) 

where  s(t)  is  a  function  of  the  intensity  modulation  impressed  upon  the 
transmitted  laser  beam.  The  simultaneous  estimators  t  and  p  may  be 
found  by  substituting  Eqs  (2)  and  (5)  into  the  simultaneous  equations 


3  £n  [f({t^},  N^It,p)] 


3  An[f ({t^},N^lT,p)] 
3^"  ; 


Unhappily,  the  simultaneous  solutions  of  Eqs  (6)  and  (7)  are 
inseparable  and  there  is  no  solution  for  t  independent  of  p  and 
visa-versa  [3,pp.l4-15] . 
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However,  substituting  Eqs  (2)  and  (5)  into  Eq' (3)  yields 


r 


Maximize  hln[f ({t^},N^|T,p)]  = 
on  T , p  ^ 


m 

E  iln[ps(t.-T)  +  X  ] 
i=l  ^  " 


-  J  [ps(a-T)  +  X^]daJ 


(8) 


The  general  form  of  Eq  (8)  is  that  of  a  correlator  [3,pp.  15-16] 
and  may  be  realized  approximately  by  a  m  x  n  bank  of  correlators 
(Fig  1).  The  "biggest  picker"  picks  the  output  of  the  correlator 
with  the  largest  amplitude,  and  the  ,p .  and  t.  that  corresponds  to  that 
correlator  are  the  estimates  p  and  t  respectively.  The  correlator 
realization  shown  in  Fig.  1  requires  in  x  n  correlators  and  is  not 
practical  to  implement . 

It  can  be  shown  that  in  the  case  of  a  background  noise  limited 
detector,  Eq  (8)  may  be  approximated  by 


Maximize 
on  p,T 


r 

j  i(t)s(t)dt 

_''o  "  n 


ps(a-T)da 


(9) 


where  q  is  the  charge  of  an  electron  and  i(t)  is  the  output  current  of 
the  photodetector.  The  first  term  in  Eq  (9)  describes  a  correlator.  If 
the  observation  interval  [0,T]  is  much  greater  than  the  largest  ex¬ 
pected  delay  t,  i.e.,  if  end  effects  can  be  ignored,  then  the  last  term 
in  the  brackets  may  be  ignored.  Thus  in  that  case  the  detector  current 
i(t)  is  correlated  with  a  linear  function  of  s(t),  as  in  familiar  radar 
applications.  It  can  be  shown  that  dorrelators  and  matched  filters  are 
equivalent  methods  for  synthesizing  the  same  ML  estimator  [4:pp.  315- 
317].  An  alternate  implementation,  typically  used  with  sinusoidal 
signals,  is  accomplished  with  a  closed-loop  estimator  known  as  a  phase- 
locked-loop  (PLL) . 


A  practical  (suboptimum)  estimator  realization  is  shown  in  Fig.  2. 
The  delay  between  the  transmission  of  the  signal  and  the  occurrence  of 
the  peak  is  equal  to  t.  The  amplitude  of  the  peak  detector  output, 
given  the  delay  estimate  t,  is  proportional  to  the  reflectance  estimate 

P  (this  is  a  suboptimum  estimate  of  p  because  of  the  requirement  for 
knowing  t) . 

Note  that,  in  general,  the  detector  output  is  not  correlated  with 
s(t),  as  would  be  done  in  a  conventional  radar,  but  Instead  is  correlated 
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with  the  weighted  natural  log  of  ps(t)  +  Thus  it  must  be  stressed 

that  classical  radar  techniques  do  not  always  apply  to  optical  systems 
involving  direct  (power)  detection.  Only  in  the  background  noise  limited 
case,  or  when  the  signal  rate  function  is  so  large  that  the  Poisson 
statistics  become  Gaussian,  may  the  familiar  radar  ML  estimator,  a 
matched  filter,  be  used  optimumly. 


Performance  Evaluation.  It  was  stated  earlier  that  it  is  desirable 
that  var  [x-t]  and  var  [p-p]  be  small.  These  variances  are  very  dif¬ 
ficult  to  compute  in  general;  therefore,  a  lower  bound  for  the  variances 
will  be  computed  using  the  Cramer -Rao  (CR)  lower  bound.  The  CR  lower 
bound  for  the  variances  of  the  estimator  errors  for  t  and  p  may  be 
shown  to  be, 


2 

a  > 

T 


-,E 


_  2 

3£n.f({t^>,N 


T 


2  2 

p  S  (t-x)  dt 
ps(t-x)+X 

n 


2  r  3inf({t,},N„Ix,p)’ 

»  •  ®  ap2 

•*  * 


(10) 


'(t-x) 


ps(t-x)  +  X 

n 


dt 


-1 


(11) 


[3:pp.  19-20].  Now  let  the  energy,  E  ,  contained  in  the  signal  s(t)  and 

o  S  . 

the  mean  square  bandwidth,  of  s(t)  be  defined  as 

s 

(12) 


-/ 

V  r\ 


(t)dt 


b2 

s 


J*  a)2s(a))da) 

pm 

I  S(a))da) 

v,_oo 


(13) 
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where  S  (a>)  is  the  power  spectral  density  of  s  (t) .  For  the  background 
noise  limited  case,  i.e.,  X  »ps(t),  substituting  Eqs  (12)  and  (13)  into 
Eqs  (10)  and  (11)  and  simplifying  yields 


T  p?E 


n 


S-BPF®s 


(lA) 


n 


> 

P  “  E 


(15) 


where  the  subscript  BPF  denotes  "bandpass  filtered,' 
signal  component  has  been  removed. 


i.e. 


the  DC 


Eqs  (lA)  and  (15)  are  the  primary  equations  used  to  compare  the 
various  spatial  and  temporal  coding  combinations  in  the  next  section  of 
this  paper.  For  each  different  spatial-temporal ■ coding  scheme,  the 
apparent  signal  ener^  and  mean  square  bandwidth  is  calculated  and  sub¬ 
stituted  into  Eqs  (lA)  and  (15).  The  exception  is, again, purely  sinusoidal 
modulation.  In  that  case,  s(t)  =  [1  +  Sin(2irf  t)]  is  substituted 

,2  “  9  ■  .  9 

directly  into  Eqs  (10)  and  (11)  to  obtain  bounds  on  and  a^.  The  modu¬ 
lation  format  with  the  lowest  and  s^  is  designated  the  optimum  modu¬ 
lation  from  an  estimation  accuracy  point  of  view.  It  should  be  emphasized 
that  the  CR  lower  bound  is  an  optimistic  perfoirmance  measure  and  may  not 
be  achievable  in  practice.  It  indicates  that,  regardless  of  how  good  the 
signal  processor,  the  smallest  error  possible  occurs  when  EqS  (lA)  and 
(15)  are  achieved  with  equality.  Such  a  processor  is  called  an  efficient 
estimator.  Any  suboptimum  process  will  result  in  error  variances  larger 
than  the  lowest  limits  given  by  Eqs  (lA)  and  (15). 

Another  important  performance  measure  is  maximum  unambiguous  range. 

As  stated  earlier,  the  range  is  derived  from  the  delay  estimate,  the 
relationship  between  range  and  delay  is  expressed  by  r  =  ct/2  where  c 
is  the  speed  of  light.  However,  since  the  transmitted  signal  is  periodic 
with  period  T,  a  signal  return  for  s(t-T)  looks  just  like  a  return  for 
s(t-(T+nT)  where  n  is  an  integer.  This  Implies  r  =  c(t  modulo  T)/2. 

Thus  the  maximum  unambiguous  range  that  may  be  measured  is  given  by 
^amb  ~  cT/2.  It  is  highly  desirable  that  is  greater  than  the 

largest  terrain  deviation  or  target  height  to  be  measured;  this  criteria 
as  well  as  estimation  accuracy  is  used  to  choose  desirable  modulation 
formats. 
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Spatial  and  Temjporal  Coding  Analysis 

The  performance  of  the  various  spatial*'temporal  coding  combinations 
are  expressed  in  terms  of  common  constraints.  It  is  assumed  that,  unless 
so  stated,  the  transmitter -and  receiver  FOVs  are  matched,  i.e.,  subtend 
the  same  ground  area  (a  pixel)  and  are  spatially  orthogonal  from  channel 
to  channel.  The  number  of  pixels  per  scan  line  and  the  pixel  size  are 
fixed.  The  beam  is  modeled  as  moving  discretely  from  plxel-to-plxel  with 
no  spatial  filtering.  All  dwell  times  are  expressed  in  terms  of  that  for 
a  single  source  laser' line  scan  system,  t^.  The  total  number  of  sources 
will  be  N,  regardless  of  the  spatial  configuration.  Noise  considerations 
are  limited  to  background  noise  and  noise  associated  with  the  direct  de^ 
tection  process.  Ideal  optics  and  atmosphere  are  assumed.  And,  of  most 
importance,  the  power  incident  upon  the  detector,  P^,  is  assumed  to  be  a 
constant,  regardless  of  the  Imposed  intensity  modulation.  Finally,  the 

signal  processor  is  always  assumed  to  be  the  appropriate  ML  estimator 
as  discussed  earlier.. 

Many  combinations  of  temporal  and  spatial  codes  were  analyzed  and  are 
reported  in  [3].  The  major  temporal  intensity  codes  included  continuous 
wave  (CW)  modulation,  sinusoidal  modulation,  pulse  code  modulation  (PCM), 
and  On-Off-Keyed  (OOK)  modulation  (Pig.  3).  As  shown  in  Fig.  3,  PCM 
modulation  is  merely  intensity  modulation  with  a  particular  pattern  of 
on-off  pulses,  each  of  width  t  .  The  pulse  pattern  was  chosen  for  its 
desirable  autocorrelation  and  crosscorrelation  properties;  one  such  well- 
known  pattern  is  the  pseudo-noise  (PN)  sequence.  Finally,  OOK  modulation 
is  the  combination  of  sinusoidal  and  PCM  modulations,  i.e.,  the  envelope 
of  the  sinusoid  is  a  PCM  pulse  sequence.  The  advantages  of  each  is  dis¬ 
cussed  in  a  later  paragraph. 

The  spatial  codes  or  scan  configurations  included  staring  (non-scanning) 
systems,  serially  scanned  multiple  beams,  multiple  beams  scanned  in  par¬ 
allel,  overlapping  scanning  beams,  and  a  staring  receiver/scanning 
multiple  sources  system.  All  of  the  logical  spatial-temporal  combinations 
possible  with  the  above  codes  are  too  numerous  to  be  reported  in  this 
paper.  For  the  sake  of  brevity,  only  the  four  most  promising  combinations 
are  discussed  here.  ^ 

The  four  most  promising  spatial-temporal  coding  combinations  are: 

(1)  Parallel  Scan,  Sinusoidal  Modulation,  (2)  Parallel  Scan,  OOK 
Sinusoidal  Modulation,  (3)  Serial  Scan,  PCM,  and  (4)  Overlapping  Beams, 

PCM.  The  first  two  spatial-temporal  coding  configurations  employ  multiple 
sources  scanned  in  parallel  (Fig  4).  By  scanning  N  lines  at  once,  the 
scan  rate  may  be  reduced  by  a  factor  of  N~^,  thereby  Increasing  the  dwell 
time,  and  consequently,  the  reflected  energy  per  pixel  by  a  factor  of  N. 

The  third  coding  configuration  employs  serially  scanned  multiple  sources 
(Fig  5) .  By  combining  the  output  of  the  N  matched  filters  via  a  delay 
line,  the  effect  is  the  same  as  looking  at  a  given  pixel  N  times,  thereby 
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increasing  the  effective  energy  reflected  by  that  pixel  by  a  factor  of  N. 
An  advantage  of  the  serial  scan  scheme  is  that  any  differences  in  trans¬ 
mitter-receiver  pair  responsivities  are  averaged  out  by  the  integrating 
process.  The  last  coding  configuration  utilizes  N  scanning,  broadly 
overlapping  transmit  beams  partitioned  by  N  scanning  receiver  (Fig  6). 

Once  again  the  effect  is  to  raise  to  energy  reflected  per  pixel  by  N. 

This  scheme  also  has  the  advantage  of  averaging  out  any  differences  in 
the  N  transmitter  powers  as  well  as  alleviating  the  problem  of  collimating 
the  relatively  broadly  diverging  GaAs  laser  beams. 

The  root  mean  square  (rms)  errors  of  the  above  systems  are  compared 
in  Table  1.  For  each  coding  scheme,  E  ,  E  _  j  and  were  Calculated  in 

S'  .  S 

terms  of  the  common  system  parameters.  Then,  in  order  to  facilitate  the 
system  comparisons,  the  rms  range  and  reflectance  performance  for  each 
system  is  expressed  in  terms  of  that  of  the  system  with  the  smallest 
mean  square  error.  The  smallest  mean  square  range  error  occurred  for 
the  Parallel  Scan,  Sinusoidal  Modulation  case  and  is  given  by 


2Tr^f  ■^p^S^p^T , 
.  m  d  d 


(16) 


where  c  is  the  speed  of  light,  f^  is  the  sinusoidal  Intensity  modulation 
frequency,  and  is  the  signal  power  Incident  upon  the  detector.  The 
smallest  mean  square  reflectance  error  occurred  for  the  two  schemes  em¬ 
ploying  PCM,  and  is  given  by 


a2  >  i 


p-PCM  -  N  '  "FpitT 
da 


There  are  several  important  observations  which  should  be  made  about 
the  results  tabulated  in  Table  1.  First,  note  that  in  order  for  the  rms 
range  error  for  the  PCM  systems  to  be  lower  than  that  for  the  sinusoidally 
modulated  system,  the  PCM  pulse  width  t^  (Fig  3)  must  be  shorter  than 

the  period  of  f^^j.  Thus  for  a  fixed  system  bandwidth  (determined  either 
by  detector  or  laser  modulation  limitations) ,  the  sinusoidal  system 
achieves  a  finer  range  accuracy.  However,  the  rms  range  error  and 
maximum  unambiguous  range  are  both  Inversely  proportional  to  the 
Intensity  modulation  frequency  f^.  Thus  range  error  performance  and 

unambiguous  range  are  direct  trade-offs. 

A  compromise  exists  in  the  OOK  technique.  The  same  sinusoidal  in¬ 
tensity  modulation  as  used  above  is  coherently  keyed  on  and  off  by  a 
PCM  code.  The  sinusoidal  modulation  Inside  of  each  burst  or  envelope 
may  be  used  for  fine  grain  range  measurement,  and  the  PCM  envelope  may 
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be  used  for  coarse  long  range  measurement.  A  small  price  is  paid  in 
range  error  performance,  i.e.,  is  increased  by  /2,  but  the  unambiguous, 
range  capability  is  significantly  increased.  There  is  also  a  similar 
degradation  of  Op,  but  the  system  will  probably  he  capable  of  Imaging 
well  after  the  signal-to-noise  ratio  is  too  low  to  make  range  estimates, 
so  this  is  not  thought  to  be  a  serious  degradation. 

The  OOK  modulation  (and  PCM)  has  one  other  significant  advantage  over 
the  sinusoidal  modulation.  By  choosing  the  proper  PCM  modulation  for  each 
source,  all  adjacent  channels  may  be  made  temporally  orthogonal  as  well  as 
spatially  orthogonal,  thereby  assuring  that  the  pixel-to-pixel  cross-talk 
will  be  an  absolute  minimum. 

Conclusions 

It  has  been  shown  that,  based  upon  the  quantum  statistics  of  a  back¬ 
ground  noise  limited  direct  detector,  the  optimum  estimator  for  delay  (or 
range),  is  a  correlator,  and  a  suboptimum  estimator  for  reflectance  is  an 
amplitude  measurement,  given  the  delay.  The  best  possible  range  performance 
is  achieved  using  multiple  sources  scanned  in  parallel  and  sinusoidally 
intensity  modulated.  However,  if  a  large  unambiguous  range  capability 
is  important,  a  good  compromise  is  a  system  using  multiple  sources  scanned 
in  parallel  and  intensity  modulated  by  a  PCM  On-Off-Keyed  sinusoid.  This 
modulation  format  also  gives  an  additional  advantage  of  using  PCM  pulse 
sequences  for  adjacent  channels  so  that  they  are  temporally  as  well  as. 
spatially  orthogonal. 
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Implementation  of  ML  Estimators 


Figure  2  Suboptimum  ML  Estimators 
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Sinusoidal  Moduation 
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Vulos  Coda  Modulation 


On>Ofr*Koyod  Modulation 


Figure  3  Primary  Intensity  Modu¬ 
lation  Codes 
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Figure  5  N  Sources  in  Series 


Laser  Beams 
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Table  I 

Spatial-Temporal  Systems  Characteristics 

Range 


System 

a  . 
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P 

Ambiguity 

Interval 

,  Optics 
Complexity 

Serial  Scan,  PCM 

2  m*'o^r-s 

%-PCM 

Overlapping  Beams,  PCM 

iVo^^r-S 

‘^p-PCM 

CTd/2 

Least 

Parallel  Scan,  Single 
Sinusoidal  Modulation 

‘^r-S 

^•1%-PCM 

Severe 

Patallel  Scan,  OOK 
Sinusoidal  Modulation 

/2  a 

r-s 

’■•“'•p-PCM 

CTd/2 
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